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W81XWH-1 2-1 -0604  (OC11Q435)  Annual  Report,  Oct  2013 

Ph  Brad  H.  Nelson,  Ph.D. 

CctRlsr  Rob  Holt,  Ph.D.,  John  Webb  Ph.D.,  Peter  Watson,  M.D. 

Title  of  Project  Deciphering  the  Adaptive  Immune  Response  to  Ovarian  Cancer 

INTRODUCTION: 

Tumor-infiltrating  CD8+  T  cells  are  strongly  associated  with  increased  survival  in  ovarian 
cancer.  However,  they  do  not  work  in  isolation.  We  discovered  that  two  other  types  of  immune 
cell  play  an  important  supportive  role:  B  cells  and  helper  T  cells  (specifically,  helper  T  cells  that 
express  a  protein  called  FoxP3).  We  made  this  discovery  by  performing  a  systematic  analysis  of 
immune  cells  in  ovarian  cancer.  We  found  that  killer  T  cells  are  often  found  in  small  clusters 
together  with  B  cells  and  helper  T  cells.  Importantly,  we  found  that  patients  whose  tumors  have 
these  combinations  of  immune  cells  have  better  survival  rates  than  patients  whose  tumors 
contain  killer  T  cells  alone.  This  tells  us  that  T  cells  and  B  cells  work  together  to  attack  tumors. 
These  findings  have  powerful  clinical  implications:  to  enhance  the  immune  response  to  ovarian 
cancer,  we  need  to  enhance  the  activity  of  all  three  types  of  immune  cell,  rather  than  killer  T 
cells  alone. 


To  explain  these  observations,  we  hypothesized  that  B  cells  serve  as  “organizers”  that 
help  to  draw  T  cells  into  the  tumor.  In  addition,  B  cells  might  present  tumor  proteins  to  the  T 
cells  to  facilitate  tumor  recognition.  We  further  hypothesized  that  FoxP3  helper  T  cells  might 
produce  cytokines  that  help  to  excite  the  killer  T  cells.  To  test  these  hypotheses,  we  proposed  to 
determine  which  tumor  proteins  (antigens)  are  recognized  by  B  cells  and  FoxP3  helper  T  cells  in 
ovarian  cancer.  By  identifying  these  antigens,  we  will  be  able  to  create  new  molecular  tools  to 
elucidate  how  the  immune  system  recognizes  and  attacks  ovarian  cancer.  The  study  has  four 
tasks: 


Task  1.  To  identify  tumor  antigens  recognized  by  CD20+  TIL. 

Task  2.  To  identify  tumor  antigens  recognized  by  CD4+FoxP3+  TIL. 

Task  3.  To  determine  whether  tumor-infiltrating  B  cells  and  T  cells  recognize  the  same 
antigens. 

Task  4.  To  assess  the  functional  phenotype  of  antigen-specific  CD4+FoxP3+  TIL. 

Significance:  The  immune  system  has  a  profound  influence  on  survival  from  ovarian 
cancer.  With  better  understanding  of  the  immune  response,  it  will  be  possible  to  design  new 
treatments  such  as  vaccines  that  enhance  tumor  immunity  and  increase  patient  survival.  We 
envision  our  work  will  lead  to  a  major  re-think  about  cancer  vaccines:  instead  of  simply  trying  to 
activate  killer  T  cells,  we  also  need  to  find  effective  ways  to  activate  their  team  mates,  the  B 
cells  and  FoxP3+  helper  T  cells. 

BODY: 

Task  1.  To  identify  tumor  antigens  recognized  by  CD20+  TIL. 

In  this  task,  we  proposed  to  identify  the  antigens  recognized  by  the  3  most  abundant 
CD20+  TIL  clones  from  each  of  3  ovarian  cancer  patients.  To  accomplish  this,  we  proposed  to 
clone  immunoglobulin  G  (IgG)  molecules  from  individual  CD20+  TIL.  These  will  be  used  to 
identify  the  corresponding  antigens  using  three  different  approaches:  candidate  antigen  assays, 
cDNA  library  screening,  and  mass  spectrometry. 
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Progress  to  date: 

Generating  IgG  anc[  TCR  /profiles  from  ovarian  tumor  specimens 

To  identify  the  3  most  abundant  IgG  (for  Task  1 )  and  TCR  clones  (for  Task  2)  in  each 
patient’s  tumor,  we  are  using  deep  sequencing  methods  developed  in  the  Holt  lab.  Our 
approach  for  deep  sequence  analysis  of  T  cell  receptor  (TCR)  beta  chain  diversity  has  been 
previously  described  (1 , 2),  and  we  have  recently  made  adaptations  to  enable  B  cell  receptor 
(BCR)  heavy  chain  profiling  and  sample  multiplexing,  as  described  below.  For  the  present 
investigation,  we  isolated  total  RNA  from  bulk  ascites  from  the  three  HGSC  ovarian  cancer 
study  subjects  at  three  different  timepoints  (the  time  of  primary  disease,  at  first  recurrence,  and 
at  second  recurrence).  First  strand  cDNA  was  synthesized  from  total  RNA  using  either  TCR 
beta  chain  or  BCR  heavy  chain  gene-specific  primers,  and  this  cDNA  was  then  used  as 
template  to  amplify  the  CDR3  sequences  from  each  receptor  population  by  PCR.  The  CDR3 
region  is  the  site  of  VDJ  recombination  and  is  the  most  highly  variable  and  thus  most 
informative  region  for  the  purpose  of  profiling  repertoire  diversity.  The  PCR  primers  were  tagged 
with  unique,  6  base  pair,  error  tolerant  barcode  sequences.  A  different  barcode  was  used  for 
each  individual  sample.  The  CDR3  amplicons  were  then  subjected  to  several  additional  rounds 
of  nested  PCR  using  primers  tailed  with  lllumina  adapter  sequences  necessary  for  annealing  to 
the  lllumina  flow  cell  surface  for  massively  parallel  sequencing.  The  samples  were  pooled  and 
sequenced  using  a  single  flow  cell  and  the  lllumina  miSeq  platform.  Sequence  data  was  then 
de-convoluted  based  on  the  barcode  identifiers,  and  the  TCR  and  BCR  CDR3  sequences  were 
filtered  to  remove  low  quality  reads,  then  clustered  using  standard  bioinformatic  methods. 
Clustering  was  done  to  determine  the  number  of  distinct  CDR3  sequences  present  in  each 
sample,  which  reflects  the  diversity  of  distinct  T-cell  and  B-cell  clonotypes  with  the  original 
ascites  specimens.  Prior  to  reverse  transcription  and  PCR  each  sample  was  spiked  with  RNA 
from  a  single,  known  clonotype  in  order  to  monitor  PCR  and  sequence  error  rates  and  the 
efficiency  of  sequence  recovery. 

So  far,  we  have  filtered  and  analyzed  the  sequence  data  derived  from  specimens  taken 
from  each  of  the  three  timepoints  for  one  of  the  three  study  subjects.  We  obtained  over  30,000 
total  TCR  beta  chain  CDR3  sequences  from  each  of  these  three  samples,  and  over  15,000  total 
BCR  heavy  chain  CDR3  sequences  from  each  of  these  same  three  samples.  These  collapsed, 
upon  clustering,  into  a  much  smaller  number  of  distinct  sequences.  Specifically  there  were 
1 ,598  distinct  TCR  sequences  and  688  distinct  BCR  sequences  observed  per  sample. 

However,  as  expected,  the  majority  of  these  unique  sequences  were  “singletons”,  representing 
rare  clonotypes,  with  the  bulk  of  the  accumulated  sequence  data  from  each  sample  represented 
approximately  30  clonotypes  per  sample.  Thus  the  ascites  repertories  are  polyclonal,  but  not  to 
the  same  extent  as  peripheral  blood,  where  previously  we  have  been  able  to  detect  over  a 
million  distinct  sequences  from  a  single  sample  by  deep  sequencing.  Interestingly,  18  of  the 
BCR  clonotypes  were  detectable  at  all  three  timepoints,  albeit  with  varying  abundance,  but  the 
TCR  clonotypes  were  much  more  sample  specific,  with  only  one  TCR  clonotype  being  present 
in  all  three  samples. 

Amplification  of  matched  IgG  heavy  anc[  light  chains  by  single-cell  RT-PCR 

Single  cell  sorting  and  PCR  amplification  of  matched  immunoglobulin  (Ig)  heavy  and 
light  chains  from  tumor-infiltrating  B  lymphocytes  is  being  accomplished  via  the  methods  of  Tiller 
et  al.  (3).  Single  B  cells  are  isolated  based  on  their  surface  expression  of  CD19,  CD20  and  IgG 
using  a  BD  Influx  cell  sorter.  Individual  cells  are  sorted  directly  into  single  wells  of  96-well  PCR 
plates.  Reverse  transcription  is  accomplished  directly  within  the  sorting  plate,  which  yields 
cDNA  corresponding  to  the  original  single  cell.  This  cDNA  is  then  split  equally  between  four 
separate  PCR  reactions,  one  each  for  the  variable  portions  of  Ig  heavy,  kappa,  and  lambda,  and 
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one  for  the  control  housekeeping  gene  GAPDH.  The  variable  portions  of  the  immunoglobulin 
genes  are  then  amplified  by  sequential  nested  multiplex  PCR  using  primer  sets  known  to 
amplify  all  human  variable  and  joining  gene  segments  and  containing  restriction  enzyme  cut 
sites  to  facilitate  downstream  molecular  cloning.  To  date  we  have  successfully  amplified 
matched  pairs  of  Ig  genes  from  B  cells  isolated  from  the  tumors  of  three  high-grade  serous 
ovarian  cancer  patients  (Figure  1).  Roughly  70%  of  all  sorting  wells  yield  amplification  of  at  least 
one  B  cell-specific  gene  product,  while  about  10-25%  yield  matched  pairs  of  heavy  and  light 
chains.  In  total,  we  have  successfully  amplified  28  pairs  of  Ig  genes  from  the  3  patients. 

Cloning  erf  CDR3  regions  into  IgG  expression  vectors 

Molecular  cloning  of  the  PCR  amplified  variable  portions  of  matched  heavy  and  light 
chain  genes  derived  from  single  sorted  tumor  infiltrating  B  cells  is  accomplished  by  standard 
techniques.  All  pairs  that  have  been  amplified  have  subsequently  been  cloned  into  appropriate 
expression  vectors  containing  the  signal  sequence  and  constant  portions  of  the  corresponding 
germline  immunoglobulin  chain  (heavy,  kappa,  lambda).  Once  inserted  into  these  expression 
vectors,  the  DNA  sequence  of  each  immunoglobulin  chain  is  determined  by  Sanger  sequencing. 
Sequences  for  each  of  the  pairs  of  heavy  and  light  chains  have  been  obtained  (Figure  2).  We 
are  currently  evaluating  the  prevalence  of  each  immunoglobulin  heavy  sequence  by  sequencing 
samples  prepared  by  conventional  RT-PCR  of  bulk  tumor  preparations. 

Preparation  of  recombinant  IgG  [ rlgG ) 

To  produce  recombinant  IgG,  two  expression  vectors  encoding  cloned  matched  heavy 
and  light  chain  variable  portions  are  co-transfected  into  log-phase  293T  cells.  This  transient 
transfection  is  accomplished  by  a  calcium  phosphate  precipitation  method.  After  seven  days  in 
culture,  the  supernatant  is  harvested  and  assessed  for  the  presence  of  rlgG  by  human  IgG- 
specific  western  blot.  Fully  assembled  IgG  molecules  (~150  kDa)  are  observed  in  the 
supernatant  of  80-90%  of  culture  transfected  with  different  pairs  of  expression  vectors  (Figure  3) 
under  non-reducing  conditions.  To  isolate  rlgG  from  the  culture  supernatant,  it  is  concentrated 
and  run  over  protein  G  sepharose  columns.  Captured  IgG  molecules  are  eluted  at  low  pH  and 
the  resulting  eluate  is  neutralized  and  desalted.  Further  concentration  and  assessment  of  purity 
is  accomplished  by  BCA  assay  and  additional  western  blot  analysis.  To  date,  1 1  recombinant 
antibodies  have  successfully  been  expressed  in  vitro  from  14  different  matched  Ig  pairs.  Based 
on  this,  we  expect  to  have  an  80%  success  rate  in  the  production  of  rlgG  from  the  remaining 
cloned  Ig  pairs. 

Task  2.  To  identify  tumor  antigens  recognized  by  CD4+FoxP3+  TIL. 

We  proposed  to  clone  TCR  molecules  from  CD4+FoxP3+  tumor-infiltrating  T  cells  from 
the  3  ovarian  cancer  patients  described  above.  Cloned  TCR’s  will  be  expressed  in  a  hybridoma 
cell  line,  which  in  turn  will  be  used  to  identify  the  underlying  antigens  using  a  candidate 
approach  and/or  cDNA  library  screening. 

Progress  to  date: 

Amplification  of  paired  TCR-alpha  anc[  - beta  CDR3  regions  by  single-cell  RT-PCR 

We  are  developing  single-cell  analytical  methods  for  profiling  CDR3  sequences  from 
paired  alpha  and  beta  subunits.  To  date,  large-scale  T-cell  repertoire  analysis  has  been  limited 
to  interrogation  of  a  single  TCR  subunit,  usually  the  beta  chain,  per  sequencing  run  (as 
described  for  Task  1 ).  Recently,  however,  the  potential  for  pairwise  lllumina  sequencing  of 
a|3TCRs  has  been  demonstrated,  by  means  of  alpha  and  beta  chain  fusion  PCR  executed 
within  single  cell  containing  droplets  of  water  in  oil  emulsions  (4).  Using  this  approach,  hundreds 
of  a|3TCR  sequences  could  be  identified  from  starting  populations  of  greater  than  1  million  total 
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cells.  Likewise,  it  was  recently  demonstrated  that  thousands  of  immunoglobulin  heavy  and  light 
chain  pairs  can  be  obtained  by  bead  capture  of  single  B-cell  mRNA  followed  by  linkage  PCR  in 
single  bead-containing  emulsion  droplets  (5).  While  the  yield  of  each  of  these  approaches  is 
modest,  they  represent  important  advances  towards  the  goal  of  deep,  cheap  and  fast  profiling  of 
dimeric  antigen  receptors.  We  have  been  pursuing  pairwise  sequencing  of  a|3TCRs  and  have 
established  a  unique  method  for  fusion  of  a|3TCR  cDNAs  during  reverse  transcription.  We  have 
demonstrated  this  approach  using  bulk  RNA  from  Jurkat  cells,  which  are  a  T  cell  line  expressing 
a  defined  a|3TCR  (Figure  4).  We  are  working  to  adapt  this  system  to  a  single  cell  microfluidics 
platform  in  order  to  establish  a  robust  single  cell  assay  amenable  to  deep  sequencing  on  the 
lllumina  platform.  In  future,  this  assay  will  allow  us  to  move  from  deep  TCR  sequence  profiling 
data  directly  to  receptor  reconstitution  (using  TCR  expression  vectors  we  now  have  in  hand, 
generously  provided  by  our  colleague  Klaus  Dornmair,  Univ.  of  Munich)  and  then  antigen 
screening.  Thus,  we  are  making  excellent  progress  in  terms  of  developing  the  necessary 
technology  to  complete  Task  2. 

Task  3.  To  determine  whether  tumor-infiltrating  B  cells  anc[  T  cells  recognize  the  same 
antigens. 

Once  we  have  identified  cognate  antigens  for  the  predominant  CD8+  (from  the  original 
IDEA  proposal),  CD20+  (Task  1)  and  CD4+FoxP3+  (Task  2)  TIL,  we  proposed  to  assess  the 
extent  to  which  these  antigen  sets  overlap.  If  the  antigens  recognized  by  CD20+  TIL  are 
also  recognized  by  CD8+  or  CD4+  TIL,  this  would  support  our  hypothesis  that  CD20+  TIL 
can  serve  as  APC  in  the  tumor  environment. 

Progress  to  date: 

This  task  is  not  scheduled  to  start  until  Year  2. 

Task  4.  To  assess  the  functional  phenotype  of  antigen-specific  CD4+FoxP3+  TiL. 

CD4+FoxP3+  TIL  show  great  functional  heterogeneity  in  EOC,  which  makes  it  difficult  to 
draw  definitive  conclusions  about  their  role  in  the  tumor  environment.  With  knowledge  of  their 
cognate  antigens,  we  will  be  able  to  clarify  this  issue  by  assessing  the  functional  phenotype  of 
individual  CD4+FoxP3+  T  cell  clones  as  opposed  to  bulk  cell  preparations.  We  proposed  to  do 
this  by  constructing  MHC  class  II  tetramers,  which  bind  specifically  to  CD4+  T  cells  expressing 
TCRs  relevant  to  a  particular  antigen. 

Progress  to  date: 

This  task  is  not  scheduled  to  start  until  Year  2. 


KEY  RESEARCH  ACCOMPLISHMENTS: 

1 .  We  published  a  key  manuscript  describing  the  results  of  the  original  IDEA  award:  Castellarin 
M,  Milne  K,  Zeng  T,  Tse  K,  Mayo  M,  Zhao  Y,  Webb  JR,  Watson  PH,  Nelson  BH,  Holt  RA. 
Clonal  evolution  of  high-grade  serous  ovarian  carcinoma  from  primary  to  recurrent  disease. 

J  Pathol.  2013  Mar;229(4):51 5-24.  doi:10.1002/path.4105.  Epub  2012  Nov  29.  PubMed 
PMID:  22996961.  (Appendix  A) 

•  A  second  key  manuscript  from  the  original  IDEA  award  is  under  revision  at  Clinical  Cancer 
Research  (with  an  enthusiastic  invitation  to  re-submit):  Wick,  D.,  Webb,  J.R.,  Nielsen,  J.S., 
Martin,  S.,  Kroeger,  D.R.,  Milne,  K.,  Castellarin,  M.,  Twumasi-Boateng,  K.,  Watson,  P.H., 
Holt,  R.H.,  Nelson,  B.H.  2013.  Surveillance  of  the  tumor  mutanome  by  T  cells  during 
progression  from  primary  to  recurrent  ovarian  cancer.  Clin  Cancer  Res,  in  revision. 
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•  As  listed  below,  in  2012-2013  our  team  published  9  other  manuscripts  with  direct  relevance 
to  this  project,  and  16  with  indirect  relevance. 

•  Task  1:  lllumina-based  IgG  profiling  has  been  completed  for  3  patients  at  3  time  points. 

•  Task  1:  Robust  methods  have  been  developed  to  amplify  matched  IgG  heavy  and  light 
chains  by  single-cell  RT-PCR. 

•  Task  1 :  Robust  methods  have  been  developed  to  clone  CDR3  regions  into  IgG  expression 
vectors. 

•  Task  1 :  Robust  methods  have  been  developed  to  prepare  recombinant  IgG  (rlgG). 

•  Task  2:  lllumina-based  profiling  of  TCRs  has  been  completed  for  3  patients  at  3  time  points; 
the  next  step  is  to  do  this  for  the  CD4+CD25+  subset. 

•  Task  2:  Robust  methods  have  been  developed  to  amplify  paired  TCR-alpha  and  -beta 
CDR3  regions  by  single-cell  RT-PCR;  the  next  step  is  to  insert  these  into  expression  vectors 
and  transfect  murine  T  cell  hybridoma  cells. 


REPORTABLE  OUTCOMES: 

Manuscripts  with  direct  relevance  published  by  the  team  in  2012/2013: 

1 .  Castellarin  M,  Milne  K,  Zeng  T,  Tse  K,  Mayo  M,  Zhao  Y,  Webb  JR,  Watson  PH,  Nelson  BH, 
Holt  RA.  Clonal  evolution  of  high-grade  serous  ovarian  carcinoma  from  primary  to  recurrent 
disease.  J  Pathol.  2013  Mar;229(4):515-24.  doi:10.1002/path.4105.  Epub  2012  Nov  29. 
PubMed  PMID:  22996961 .  (Appendix  A) 

2.  Nielsen  JS,  Sahota  RA,  Milne  K,  Kost  SE,  Nesslinger  NJ,  Watson  PH,  Nelson  BH.  CD20+ 
tumor-infiltrating  lymphocytes  have  an  atypical  CD27-  memory  phenotype  and  together  with 
CD8+  T  cells  promote  favorable  prognosis  in  ovarian  cancer.  Clin  Cancer  Res.  2012  Jun 

1 5;  1 8(1 2):3281  -92.  doi:  10.1 158/1 078-0432.CCR-1 2-0234.  Epub  2012  May  2.  PubMed 
PMID:  22553348.  (Appendix  B) 

3.  deLeeuw  RJ,  Kost  SE,  Kakal  JA,  Nelson  BH.  The  prognostic  value  of  FoxP3+  tumor- 
infiltrating  lymphocytes  in  cancer:  a  critical  review  of  the  literature.  Clin  Cancer  Res.  2012 
Jun  1 ;  1 8(1 1  ):3022-9.  doi:  10.1 158/1 078-0432.CCR-1 1-3216.  Epub  2012  Apr  17.  Review. 
PubMed  PMID:  22510350.  (Appendix  C) 

4.  Nielsen  JS,  Nelson  BH.  Tumor-infiltrating  B  cells  and  T  cells:  Working  together  to  promote 
patient  survival.  Oncoimmunology.  2012  Dec  1 ;  1  (9):  1 623-1 625.  PubMed  PMID:  23264915; 
PubMed  Central  PMCID:  PMC3525624. 

5.  Watson  CT,  Steinberg  KM,  Huddleston  J,  Warren  RL,  Malig  M,  Schein  J,  Willsey  AJ,  Joy  JB, 
Scott  JK,  Graves  TA,  Wilson  RK,  Holt  RA,  Eichler  EE,  Breden  F.  Complete  haplotype 
sequence  of  the  human  immunoglobulin  heavy-chain  variable,  diversity,  and  joining  genes 
and  characterization  of  allelic  and  copy-number  variation.  Am  J  Hum  Genet.  2013  Apr 
4;92(4):530-46.  doi:  10.1016/j.ajhg.2013.03.004.  Epub  2013  Mar  28.  PubMed  PMID: 
23541343;  PubMed  Central  PMCID:  PMC3617388. 

6.  Warren  RL,  Choe  G,  Freeman  DJ,  Castellarin  M,  Munro  S,  Moore  R,  Holt  RA.  Derivation  of 
HLA  types  from  shotgun  sequence  datasets.  Genome  Med.  2012  Dec  10;4(12):95.  [Epub 
ahead  of  print]  PubMed  PMID:  23228053;  PubMed  Central  PMCID:  PMC3580435. 

7.  Milne  K,  Alexander  C,  Webb  JR,  Sun  W,  Dillon  K,  Kalloger  SE,  Gilks  CB,  Clarke  B,  Kobel  M, 
Nelson  BH.  Absolute  lymphocyte  count  is  associated  with  survival  in  ovarian  cancer 
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independent  of  tumor-infiltrating  lymphocytes.  J  Transl  Med.  2012  Feb  27;10:33.  doi: 
10.1186/1479-5876-10-33.  PubMed  PMID:  22369276;  PubMed  Central  PMCID: 
PMC3310776. 

8.  Nelson  BH.  Killer  T  cells  to  the  rescue  in  ovarian  cancer.  Gynecol  Oncol.  2012 
Feb;124(2):178-9.  doi:  10.1016/j.ygyno.201 1.12.434.  PubMed  PMID:  22264602. 

9.  Woodsworth  D.,  Castellarin  M.,  Holt  R.A.  2013.  Sequence  analysis  of  T-cell  repertoires  in 
health  and  disease.  Genome  Medicine,  in  press. 

10.  Webb,  J.R.,  Milne,  K.,  Watson,  P.H.,  deLeeuw,  R.D.,  and  Nelson,  B.H.  2013.  Tumor- 
infiltrating  lymphocytes  expressing  the  tissue  resident  memory  marker  CD103  are 
associated  with  increased  survival  in  high-grade  serous  ovarian  cancer.  Clin  Cancer  Res,  in 
press. 

Manuscripts  with  indirect  relevance  published  by  the  team  in  2012/2013: 

1.  Sio  A,  Chehal  MK,  Tsai  K,  Fan  X,  Roberts  ME,  Nelson  BH,  Grembecka  J,  Cierpicki  T,  Krebs 
DL,  Harder  KW.  Dysregulated  hematopoiesis  caused  by  mammary  cancer  is  associated 
with  epigenetic  changes  and  hox  gene  expression  in  hematopoietic  cells.  Cancer  Res.  2013 
Oct  1;73(19):5892-904.  doi:1 0.1 158/0008-5472.CAN-1 3-0842.  Epub  2013  Aug  1.  PubMed 
PMID:  23913828. 

2.  West  NR,  Kost  SE,  Martin  SD,  Milne  K,  Deleeuw  RJ,  Nelson  BH,  Watson  PH.  Tumour- 
infiltrating  FOXP3(+)  lymphocytes  are  associated  with  cytotoxic  immune  responses  and 
good  clinical  outcome  in  oestrogen  receptor-negative  breast  cancer.  Br  J  Cancer.  2013  Jan 
15;1 08(1  ):1 55-62.  doi:  10.1038/bjc.2012.524.  Epub  2012  Nov  20.  PubMed  PMID: 

23169287;  PubMed  Central  PMCID:  PMC3553524. 

3.  Cancer  Genome  Atlas  Research  Network.  The  Cancer  Genome  Atlas  Pan-Cancer  analysis 
project.  Nat  Genet.  2013  Sep  26;45(10):1 1 13-20.  doi:  10.1038/ng.2764.  PubMed  PMID: 
24071849. 

4.  Yan  Y,  Li  X,  Blanchard  A,  Bramwell  VH,  Pritchard  Kl,  Tu  D,  Shepherd  L,  Myal  Y,  Penner  C, 
Watson  PH,  Leygue  E,  Murphy  LC.  Expression  of  both  estrogen  receptor-beta  1  (ER-pi) 
and  its  co-regulator  steroid  receptor  RNA  activator  protein  (SRAP)  are  predictive  for  benefit 
from  tamoxifen  therapy  in  patients  with  estrogen  receptor-alpha  (ER-a)-negative  early  breast 
cancer  (EBC).  Ann  Oncol.  2013  Aug;24(8):1 986-93.  doi:  10.1093/annonc/mdt132.  Epub 

201 3  Apr  1 1 .  PubMed  PMID:  2357981 6. 

5.  Cancer  Genome  Atlas  Research  Network.  Comprehensive  molecular  characterization  of 
clear  cell  renal  cell  carcinoma.  Nature.  2013  Jul  4;499(7456):43-9.  doi: 

1 0.1 038/naturel 2222.  Epub  2013  Jun  23.  PubMed  PMID:  23792563;  PubMed  Central 
PMCID:  PMC3771322. 

6.  Cancer  Genome  Atlas  Research  Network.  Genomic  and  epigenomic  landscapes  of  adult  de 
novo  acute  myeloid  leukemia.  N  Engl  J  Med.  2013  May  30;368(22):2059-74.  doi: 

1 0.1 056/NEJMoal 301689.  Epub  2013  May  1.  Erratum  in:  N  Engl  J  Med.  2013  Jul 
4;369(1  ):98.  PubMed  PMID:  23634996;  PubMed  Central  PMCID:  PMC3767041. 

7.  Cancer  Genome  Atlas  Research  Network.  Integrated  genomic  characterization  of 
endometrial  carcinoma.  Nature.  2013  May  2;497(7447):67-73.  doi:  10.1 038/naturel 21 13. 
Erratum  in:  Nature.  2013  Aug  8;500(7461):242.  PubMed  PMID:  23636398;  PubMed  Central 
PMCID:  PMC3704730. 
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8.  West  NR,  Murray  Jl,  Watson  PH.  Oncostatin-M  promotes  phenotypic  changes  associated 
with  mesenchymal  and  stem  cell-like  differentiation  in  breast  cancer.  Oncogene.  2013  Apr 
15.  doi:  10.1038/onc.2013.105.  [Epub  ahead  of  print]  PubMed  PMID:  23584474. 

9.  Cancer  Genome  Atlas  Network.  Comprehensive  molecular  portraits  of  human  breast 
tumours.  Nature.  2012  Oct  4;490(7418):61-70.  doi:  10.1038/naturel  1412.  Epub  2012  Sep 
23.  PubMed  PMID:  23000897;  PubMed  Central  PMCID:  PMC3465532. 

10.  Cancer  Genome  Atlas  Research  Network.  Comprehensive  genomic  characterization  of 
squamous  cell  lung  cancers.  Nature.  2012  Sep  27;489(7417):519-25.  doi: 

10.1038/naturel  1404.  Epub  2012  Sep  9.  Erratum  in:  Nature.  2012  Nov  8;491(7423):288. 
Rogers,  Kristen  [corrected  to  Rodgers,  Kristen].  PubMed  PMID:  22960745;  PubMed  Central 
PMCID:  PMC34661 13. 

1 1 .  Spowart  JE,  Townsend  KN,  Huwait  H,  Eshragh  S,  West  NR,  Ries  JN,  Kalloger  S,  Anglesio 
M,  Gorski  SM,  Watson  PH,  Gilks  CB,  Huntsman  DG,  Lum  JJ.  The  Autophagy  Protein  LC3A 
Correlates  with  Hypoxia  and  is  a  Prognostic  Marker  of  Patient  Survival  in  Clear  Cell  Ovarian 
Cancer.  J  Pathol.  2012  Aug  27.  doi:  10. 1002/path. 4090.  [Epub  ahead  of  print]  PubMed 
PMID:  22926683. 

12.  Cancer  Genome  Atlas  Network.  Comprehensive  molecular  characterization  of  human  colon 
and  rectal  cancer.  Nature.  2012  Jul  18;487(7407):330-7.  doi:  10.1038/naturel  1252.  PubMed 
PMID:  22810696;  PubMed  Central  PMCID:  PMC3401966. 

13.  McAlpine  JN,  Porter  H,  Kobel  M,  Nelson  BH,  Prentice  LM,  Kalloger  SE,  Senz  J,  Milne  K, 
Ding  J,  Shah  SP,  Huntsman  DG,  Gilks  CB.  BRCA1  and  BRCA2  mutations  correlate  with 
TP53  abnormalities  and  presence  of  immune  cell  infiltrates  in  ovarian  high-grade  serous 
carcinoma.  Mod  Pathol.  2012  May;25(5):740-50.  doi:  10.1038/modpathol.2011.211.  Epub 
2012  Jan  27.  PubMed  PMID:  22282309. 

14.  West  NR,  Murphy  LC,  Watson  PH.  Oncostatin  M  suppresses  oestrogen  receptor-a 
expression  and  is  associated  with  poor  outcome  in  human  breast  cancer.  Endocr  Relat 
Cancer.  2012  Apr  1 0;  1 9(2):  1 81  -95.  doi:  1 0.1 530/ERC-1 1  -0326.  Print  2012  Apr.  PubMed 
PMID:  22267707. 

15.  Castellarin  M,  Warren  RL,  Freeman  JD,  Dreolini  L,  Krzywinski  M,  Strauss  J,  Barnes  R, 
Watson  P,  Allen-Vercoe  E,  Moore  RA,  Holt  RA.  Fusobacterium  nucleatum  infection  is 
prevalent  in  human  colorectal  carcinoma.  Genome  Res.  2012  Feb;22(2):299-306.  doi: 

10.1 101/gr.  126516.1 1 1.  Epub  2011  Oct  18.  PubMed  PMID:  22009989;  PubMed  Central 
PMCID:  PMC3266037. 

16.  Madhurantakam  C,  Duru  AD,  Sandalova  T,  Webb  JR,  Achour  A.  Inflammation-associated 
nitrotyrosination  affects  TCR  recognition  through  reduced  stability  and  alteration  of  the 
molecular  surface  of  the  MHC  complex.  PLoS  One.  2012;7(3):e32805.  doi: 

10. 1371 /journal. pone. 0032805.  Epub  2012  Mar  14.  PubMed  PMID:  22431983;  PubMed 
Central  PMCID:  PMC3303804. 
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Leveraged  funding: 

In  2012,  we  obtained  three  new  grants  that,  although  not  directly  relevant  to  ovarian  cancer,  are 
helping  us  to  build  our  team  and  infrastructure  in  the  cancer  immunotherapy  field: 

1 .  The  mutated  lymphoma  genome:  a  target  for  therapeutic  vaccination 
Source:  Canadian  Cancer  Society  Research  Institute 

Dates:  02/2013-01/2015  Term:  2  years 

PI:  Brad  Nelson 

The  major  goal  of  this  project  is  to  assess  the  immunogenicity  of  the  lymphoma  genome  and 
develop  personalized  therapeutic  vaccines. 

2.  Targeting  the  Breast  Cancer  Genome  with  Personalized  Therapeutic  Vaccines 
Source:  Canadian  Breast  Cancer  Foundation  BC/Yukon 

Dates:  08/2012-07/2015  Term:  3  years 

PI:  Brad  Nelson 

The  major  goal  of  this  project  is  to  use  genome  and  transcriptome  sequencing  data  to  design 
personalized  cancer  vaccines  that  target  mutations  and  destroy  spontaneous  mammary 
tumours. 

3.  Small-molecule  inhibitors  of  the  PD1-PDL1  interaction  for  treating  metastatic  cancer 
Source:  Canadian  Cancer  Society  Research  Institute 

Dates:  02/2013-01/2015  Term:  2  years 

PI:  Jeremy  Wulff,  Co-applicant:  Brad  Nelson 

The  major  goal  of  this  project  is  to  develop  small  molecules  that  enhance  anti-tumor  T  cell 
responses  by  inhibiting  the  PD1-PDL1  interaction. 

CONCLUSION: 

Overall,  this  study  is  progressing  on  schedule  and  on  budget,  with  no  major  deviations  from  the 
original  proposal.  We  have  developed  the  necessary  methods  to  complete  Tasks  1  and  2,  which 
will  enable  progress  to  Tasks  3  and  4  as  originally  scheduled.  We  published  10  relevant 
manuscripts  in  2012-2013,  and  16  with  indirect  relevance.  Additional  funding  has  been  received 
from  several  other  agencies,  enhancing  the  strength  of  our  cancer  immunology  research 
program. 


REFERENCES: 

1 .  Warren  RL,  Freeman  JD,  Zeng  T,  Choe  G,  Munro  S,  Moore  R,  et  al.  Exhaustive  T-cell 
repertoire  sequencing  of  human  peripheral  blood  samples  reveals  signatures  of  antigen 
selection  and  a  directly  measured  repertoire  size  of  at  least  1  million  clonotypes.  Genome 
Res.  2011;21:790-7. 

2.  Freeman  JD,  Warren  RL,  Webb  JR,  Nelson  BH,  Holt  RA.  Profiling  the  T-cell  receptor  beta- 
chain  repertoire  by  massively  parallel  sequencing.  Genome  Res.  2009;19:1817-24. 

3.  Tiller  T,  Meffre  E,  Yurasov  S,  Tsuiji  M,  Nussenzweig  MC,  Wardemann  H.  Efficient 
generation  of  monoclonal  antibodies  from  single  human  B  cells  by  single  cell  RT-PCR  and 
expression  vector  cloning.  J  Immunol  Methods.  2008;329:112-24. 

4.  Turchaninova  MA,  Britanova  OV,  Bolotin  DA,  Shugay  M,  Putintseva  EV,  Staroverov  DB,  et 
al.  Pairing  of  T-cell  receptor  chains  via  emulsion  PCR.  Eur  J  Immunol.  2013;43:2507-15. 
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5.  DeKosky  BJ,  Ippolito  GC,  Deschner  RP,  Lavinder  JJ,  Wine  Y,  Rawlings  BM,  et  al.  High- 
throughput  sequencing  of  the  paired  human  immunoglobulin  heavy  and  light  chain 
repertoire.  Nat  Biotechnol.  2013;31:166-9. 

APPENDICES: 

A.  Castellarin  M,  Milne  K,  Zeng  T,  Tse  K,  Mayo  M,  Zhao  Y,  Webb  JR,  Watson  PH,  Nelson 
BH,  Holt  RA.  Clonal  evolution  of  high-grade  serous  ovarian  carcinoma  from  primary  to 
recurrent  disease.  J  Pathol.  2013  Mar;229(4):515-24.  doi:10.1002/path.4105.  Epub  2012 
Nov  29.  PubMed  PMID:  22996961. 

B.  Nielsen  JS,  Sahota  RA,  Milne  K,  Kost  SE,  Nesslinger  NJ,  Watson  PH,  Nelson  BH. 
CD20+  tumor-infiltrating  lymphocytes  have  an  atypical  CD27-  memory  phenotype  and 
together  with  CD8+  T  cells  promote  favorable  prognosis  in  ovarian  cancer.  Clin  Cancer 
Res.  2012  Jun  1 5;1 8(1 2):3281 -92.  doi:  10.1 158/1 078-0432.CCR-1 2-0234.  Epub  2012 
May  2.  PubMed  PMID:  22553348. 

C.  deLeeuw  RJ,  Kost  SE,  Kakal  JA,  Nelson  BH.  The  prognostic  value  of  FoxP3+  tumor- 
infiltrating  lymphocytes  in  cancer:  a  critical  review  of  the  literature.  Clin  Cancer  Res. 

2012  Jun  1 ;  1 8(1 1  ):3022-9.  doi:  10.1 158/1078-0432.CCR-1 1-3216.  Epub  2012  Apr  17. 
Review.  PubMed  PMID:  22510350. 


SUPPORTING  DATA: 

None. 


FIGURES: 

See  following  pages. 
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Figure  1.  PCR-based  immunglobulin  (Ig)  gene  amplification  from  single  tumor- 
infiltrating  B  cells  A)  A  representative  agarose  gel  electrophoresis  of  PCR  products 
corresponding  to  Ig-heavy  (H),  Ig-kappa  (K),  Ig-lambda  (L)  and  GAPDH  (G).  B)  Amplified 
matched  Ig  heavy  and  light  chain  variable  portions  from  44  single  B  cells  (CD19+  CD20+ 
lgG+)  sorted  from  the  disaggregated  tumor  of  an  individual  high  grade  serous  ovarian 
cancer  patient.  Each  row  contains  products  derived  from  11  single  sorted  B  cells  and  one 
control  where  no  B  cell  was  sorted  (last  set  of  reactions  in  each  row).  As  indicated,  a  total 
of  11/44  matched  heavy  and  light  chain  pairs  were  observed  (A5-D11). 
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Figure  2.  Molecular  cloning,  sequencing,  and  alignment  of  paired  immunoglobulin 
genes.  Following  PCR  amplification,  pairs  of  immunoglobulin  genes  were  molecularly 
cloned  into  expression  vectors.  Inserted  portions  were  then  sequenced  by  Sanger 
sequencing.  Following  sequencing,  the  variable  portions  of  cloned  immunoglobulin 
sequences  were  aligned  to  the  germline  elements  via  the  IMGT  V-Quest  software. 
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Figure  3.  In  vitro  production  of  recombinant  IgG  from  immunoglobulin  genes 
derived  from  tumor-infiltrating  B  lymphocytes.  After  cloning  PCR-amplified  Ig 
variable  portions  into  expression  vectors  containing  the  constant  portion  of  IgH  or 
IgK/lgL,  recombinant  antibodies  (rlgG)  were  produced  in  293T  cells  co-transfected 
with  matched  heavy  and  light  chain-encoding  vectors.  The  presence  of  fully 
assembled  IgG  antibody  in  the  supernatant  was  confirmed  by  western  blot  under 
non-reducing  conditions  as  a  150kDa  band.  The  presence  of  IgG  heavy  chain 
(~50kDa)  was  confirmed  under  reducing  conditions. 
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Figure  4.  Pairwise  sequencing  of  TCR  a  and  p  subunits.  A.  For  cDNA  synthesis  we  use  a 
mix  of  TCRa  and  TCRb  gene-specific  3’  primers,  and  the  common  SMRT  template  switching 
method  (Clontech)  to  add  5’  cDNA  priming  sites.  In  step  1 ,  the  gene-specific  oligos  (a-1  and 
b-1)  prime  reverse  transcription,  which  adds  the  5’  anchor  sequence  provided  by  the  template 
switching  oligos  TS-1  and  TS-2.  The  tail  sequence  of  TS-1  (pink)  is  the  reverse  complement  of 
the  TS-2  tail  sequence  (blue).  These  complementary  tail  sequences  are  thereby  incorporated 
into  the  reverse  transcription  products.  In  step  2,  the  complementary  sequences  incorporated 
into  the  first  strand  cDNA  prime  bidirectonal  extension  to  yield  a  single  double  stranded  DNA 
fragment  comprised  of  the  two  TCR  subunits  in  tandem  that  serve  as  template  for  further  PCR 
amplification.  By  chance,  fusions  may  be  both  TCRa  or  both  TCRp,  but  these  do  not  amplify 
due  to  suppressive  self  annealing,  and  paried  a|3  chains  are  the  main  PCR  product.  B.  Using 
total  RNAfrom  Jurkat  cells  (a  human  CD4+  T  cell  line  that  expresses  a  distinct  a|3TCR)  we 
can  effectively  amplify  the  approximately  1 .3  Kb  abTCR  fusion  product.  C.  Verification  of  the 
Jurkat  a|3TCR  fusion  product  by  capillary  sequencing. 
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Abstract 

High-grade  serous  carcinoma  (HGSC)  is  the  most  common  and  fatal  form  of  ovarian  cancer.  While  most  tumours 
are  highly  sensitive  to  cytoreductive  surgery  and  platinum-  and  taxane-based  chemotherapy,  the  majority  of 
patients  experience  recurrence  of  treatment-resistant  tumours.  The  clonal  origin  and  mutational  adaptations 
associated  with  recurrent  disease  are  poorly  understood.  We  performed  whole  exome  sequencing  on  tumour 
cells  harvested  from  ascites  at  three  time  points  (primary,  first  recurrence,  and  second  recurrence)  for  three 
HGSC  patients  receiving  standard  treatment.  Somatic  point  mutations  and  small  insertions  and  deletions  were 
identified  by  comparison  to  constitutional  DNA.  The  clonal  structure  and  evolution  of  tumours  were  inferred 
from  patterns  of  mutant  allele  frequencies.  TP53  mutations  were  predominant  in  all  patients  at  all  time  points, 
consistent  with  the  known  founder  role  of  this  gene.  Tumours  from  all  three  patients  also  harboured  mutations 
associated  with  cell  cycle  checkpoint  function  and  Golgi  vesicle  trafficking.  There  was  convergence  of  germline 
and  somatic  variants  within  the  DNA  repair,  ECM,  cell  cycle  control,  and  Golgi  vesicle  pathways.  The  vast  majority 
of  somatic  variants  found  in  recurrent  tumours  were  present  in  primary  tumours.  Our  findings  highlight  both 
known  and  novel  pathways  that  are  commonly  mutated  in  HGSC.  Moreover,  they  provide  the  first  evidence  at 
single  nucleotide  resolution  that  recurrent  HGSC  arises  from  multiple  clones  present  in  the  primary  tumour  with 
negligible  accumulation  of  new  mutations  during  standard  treatment. 
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Introduction 

Each  year,  over  220  000  women  are  diagnosed  with 
ovarian  cancer  and  140  000  die  of  the  disease  (1]. 
Of  the  four  major  histological  subtypes  of  epithelial 
ovarian  cancer,  high-grade  serous  carcinoma  (HGSC) 
is  the  most  common  and  fatal  [2,3].  HGSC  tumours  are 
thought  to  arise  from  the  ovarian  surface  epithelium 
and  the  fimbrial  region  of  the  Fallopian  tubes  [4]. 

Our  knowledge  of  the  mutational  landscape  of 
HGSC  is  improving.  Approximately  13-15%  of  HGSC 
cases  are  associated  with  germline  mutations  in  BRCA1 
and  BRCA2  [5,6],  and  additional  susceptibility  genes 
are  known  to  alter  the  BRCA1/2  DNA  repair  path¬ 
way  (eg  RAD51C,  RAD51D,  BRIP1 ,  BARD1,  CHEK2, 
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MRE11A,  NBN,  PALB2,  and  RAD50  [7]).  More¬ 
over,  at  least  one- third  of  sporadic  HGSC  cases 
show  inactivation  of  the  BRCA1  or  BRCA2  path¬ 
ways  through  acquired  mutations  or  gene  silencing 
by  DNA  methylation  [8]  such  that  altogether,  homol¬ 
ogous  recombination-based  DNA  repair  is  defective 
in  about  half  of  HGSC  cases  [9].  In  addition,  the 
tumour  suppressor  gene  TP53  is  mutated  in  over  96% 
of  HGSCs  [9,10].  Other  commonly  mutated  genes  in 
HGSC  include  EAT3,  CSMD3 ,  NF1,  RBI ,  GABRA6, 
and  CDK12  [9].  Finally,  a  common  gene  fusion  involv¬ 
ing  ESRRA  and  Cllorf20  is  found  in  10-15%  of 
HGSCs  [11]. 

Treatment  of  HGSC  typically  comprises  platinum- 
based  chemotherapeutic  agents  including  cisplatin, 
carboplatin,  and  oxaliplatin.  These  compounds  exert 
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an  anti-tumour  effect  by  inducing  intra-strand  and 
inter-strand  crosslinks  in  genomic  DNA  and  are  cell 
cycle-non-specific  in  action.  Ovarian  cancer  cells 
with  inactivating  mutations  in  the  BRCA1  or  BRCA2 
genes  are  more  prone  to  the  DNA-damaging  effects 
of  platinum  compounds.  However,  platinum  chemore- 
sistance  can  arise  through  a  variety  of  mechanisms  in 
HGSC,  including  acquired  mutations  that  restore  the 
BRCA1  and  BRCA2  open  reading  frames  [12-16]. 
Additionally,  drug  resistance  in  relapsed  tumours  can 
result  from  genetic  alterations  that  affect  drug  efflux, 
drug  metabolism,  intracellular  signalling,  membrane 
receptors,  apoptotic  signalling,  and  interference  with 
cell  replication  [17]. 

There  is  increasing  appreciation  for  the  genetic 
heterogeneity  of  epithelial  cancers  within  individual 
patients  [18-23].  The  presence  of  multiple  tumour 
clones  in  a  patient  creates  the  opportunity  for  evolu¬ 
tionary  selection  during  treatment.  Indeed,  a  previous 
study  of  recurrent  HGSC  using  FISH,  SNP  arrays, 
and  array  CGH  concluded  that  chemotherapy-resistant 
clones  arise  from  minor  clones  present  in  primary 
tumours  [24].  However,  the  extent  to  which  recur¬ 
rent  tumours  further  evolve  by  accumulating  additional 
mutations  has  yet  to  be  systematically  addressed  in 
HGSC.  Here,  we  performed  whole  exome  sequenc¬ 
ing  of  genomic  DNA  from  tumour  cells  obtained 
from  ascites  of  three  HGSC  patients  at  multiple  time 
points  (primary  treatment,  first  recurrence,  and  sec¬ 
ond  recurrence).  We  inferred  tumour  clonality  and 
relapse-associated  mutations  by  analysing  trends  in 
mutant  allele  frequencies.  The  possible  functional  con¬ 
sequences  of  mutations  were  assessed  by  pathway  anal¬ 
ysis.  Furthermore,  to  gain  a  more  comprehensive  view 
of  HGSC  genetic  susceptibility  for  these  patients,  we 
investigated  the  coincidence  of  germline  and  somatic 
variants  within  gene  interaction  networks. 


Materials  and  methods 
Clinical  specimens 

Biospecimens  and  clinical  data  were  collected  with 
informed  patient  consent  through  a  prospective  study 
entitled  BC  Cancer  Agency’s  Tumor  Tissue  Repos¬ 
itory/Immune  Response  to  Ovarian  Cancer  (IROC), 
which  was  approved  by  the  Research  Ethics  Board  of 
the  BC  Cancer  Agency  and  the  University  of  British 
Columbia.  Inclusion  criteria  included  a  diagnosis  of 
HGSC,  standard  treatment  with  surgery  followed  by 
platinum-based  chemotherapy  (with  or  without  tax- 
anes),  and  biospecimen  availability.  Peripheral  blood 
mononuclear  cells  (PBMCs)  and  primary  ascites  spec¬ 
imens  were  collected  at  the  time  of  surgery.  Recur¬ 
rent  ascites  specimens  were  collected  during  palliative 
paracentesis.  Ascites  cells  were  pelleted  by  centrifuga¬ 
tion,  cryopreserved  immediately,  and  stored  in  liquid 
nitrogen. 

Copyright  ©  2012  Pathological  Society  of  Great  Britain  and  Ireland. 
Published  by  John  Wiley  &  Sons,  Ltd.  www.pathsoc.org.uk 


Tumour  cell  enrichment 

Bulk  ascites  cells  were  thawed  and  resuspended  in 
complete  RPMI  medium  (cRPMI;  RPMI  1640,  10% 
FBS,  HEPES,  L-glutamine,  p-mercaptoethanol)  at 
a  concentration  of  107  cells  per  80  pi.  To  remove 
haematopoietic  cells,  20  pi  of  anti-human  CD45 
microbeads  (Miltenyi,  Bergisch  Gladbach,  Germany) 
was  added  per  107  total  cells  and  the  mixture  was 
incubated  at  4°C  for  15  min.  Cells  were  washed 
twice  with  10  ml  of  cRPMI  and  then  passed  over  an 
LS  column  (Miltenyi)  that  had  been  pre-wetted  with 
0.5  ml  of  cRPMI.  CD45—  cells  passing  through  the 
column  were  collected,  pelleted,  and  snap  frozen  at 
— 80  °C  for  subsequent  DNA  preparation. 

Exome  library  construction  and  sequencing 

Genomic  DNA  (gDNA)  was  isolated,  acoustically 
sheared,  size-selected  by  polyacrylamide  gel  elec¬ 
trophoresis  (PAGE),  and  end-polished.  gDNA  frag¬ 
ments  were  ligated  with  Illumina  library  adapters  and 
amplified  with  Illumina  sequencing  primers  (Illumina, 
San  Diego,  CA,  USA).  The  library  was  then  enriched 
for  exon  sequences  using  the  Agilent  Sureselect  Human 
All  Exon  50  Mb  kit  [25]  according  to  the  manufac¬ 
turer’s  instructions  (Agilent  Technologies,  Santa  Clara, 
CA,  USA).  The  exome  fraction  was  subjected  to  mas¬ 
sively  parallel  sequencing  on  the  Illumina  platform. 
Between  50  and  90  million  paired  100  bp  raw  reads 
were  obtained  from  each  sample. 

Bioinformatic  analysis 

The  short  Illumina  sequence  reads  were  aligned  to 
the  human  genome  (hgl8)  using  the  Burrows -Wheeler 
Aligner  (BWA)  [26].  Single  nucleotide  variant  (SNV) 
detection  was  performed  using  Samtools  varfilter  [27] 
and  SNVMixl  [28]  software  to  identify  high-quality, 
novel,  protein-coding  variants  between  the  sample 
sequence  and  the  human  reference  sequences.  Inser¬ 
tions  and  deletions  were  detected  using  Samtools  varfil¬ 
ter  and  then  selected  as  protein-coding  variants  as 
delineated  in  the  Hgl8  reference  sequence.  Somatic 
mutations  were  determined  by  subtracting  variants 
found  in  matched  normal  samples.  The  false-negative 
rate  for  SNVs  was  reduced  by  selecting  those  that  were 
found  by  both  Samtools  and  SNVMixl  and  by  reduc¬ 
ing  mapping  errors  using  a  targeted  de  novo  assem¬ 
bler,  TASR  [29].  All  somatic  variants  were  manually 
curated  using  the  Integrative  Genome  Viewer  (IGV) 
[30].  Germline  SNVs  were  identified  from  the  exome 
data  and  were  considered  homozygous  if  they  had  a 
greater  than  80%  probability  of  being  homozygous  as 
determined  by  SNVMixl.  Genomic  coordinates  for  the 
germline  and  somatic  variants  were  converted  from  the 
hgl8  reference  genome  to  the  hgl9  reference  genome 
using  the  UCSC  Lift  Genome  Annotations  tool  [31]. 

Biological  pathway  analysis  was  performed  using 
both  the  GeneMania  plugin  [32]  in  Cytoscape 
v2.8.3  [33]  and  Ingenuity  Pathways  Analysis 
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(Ingenuity  Systems,  Inc,  Redwood  City,  CA,  USA; 
http://www.ingenuity.com)  on  genes  with  homozygous 
germline  SNVs  that  were  found  to  be  shared  in  all 
sequencing  libraries  (tumour  and  matched  normal) 
from  all  three  patients.  Gene  identifiers  that  contained 
shared,  homozygous,  germline  SNVs  were  mapped  to 
their  corresponding  objects  based  on  direct  interac¬ 
tions  and  overlaid  onto  a  global  molecular  network 
that  was  developed  from  the  Ingenuity  Knowledge 
Base  or  the  GeneMania  human  interaction  database 
version  2012-01-06,  which  contains  102  821  535 
interactions.  The  functional  interaction  analysis  of 
a  network  identified  the  biological  functions  and/or 
diseases  that  were  most  significant  to  the  molecules  in 
the  network.  The  network  molecules  associated  with 
biological  functions  and/or  diseases  in  the  Ingenuity 
Knowledge  Base  or  GeneMania  database  were  anal¬ 
ysed.  The  Fisher’s  exact  test,  as  implemented  in  the 
Ingenuity  Knowledge  Base,  was  used  to  calculate 
the  probability  that  each  biological  function  and/or 
disease  assigned  to  that  network  was  due  to  chance 
alone.  Gene  clustering  was  performed  with  JMP8 
statistical  software  [34]  using  a  hierarchical  Ward 
clustering  algorithm  based  on  somatic  SNV  allelic 
frequencies.  Gene  Ontology  (GO)  [35]  annotated 
functions  of  gene  clusters  were  predicted  based  on 
physical  interactions  and  pathway  data  sources  using 
GeneMania  [36].  GeneMania  estimates  the  false 
discovery  rate  (FDR),  the  proportion  of  positives  that 
are  expected  to  be  false  positives,  as  g -values  by 
using  the  Benjamini-Hochberg  procedure,  and  it  is 
this  number  that  we  report.  Nodes  were  partitioned 
according  to  network  modules  using  the  Reactome  FI 
Cytoscape  plugin  [37]  and  network  diagrams  were 
depicted  using  the  circular  layout  in  Cytoscape. 

Sequence  validation 

Sequence  validation  was  performed  by  PCR  amplifi¬ 
cation  of  targeted  mutant  loci  followed  by  Illumina 
sequencing  of  indexed  products  on  a  single  cell  lane  of 
100  bp  paired-end  reads  on  an  Illumina  HiSeq  platform. 
The  short  Illumina  sequence  reads  were  then  aligned  to 
the  human  genome  (hgl8)  using  the  BWA  [26].  Variant 
detection  was  performed  using  Samtools  mpileup  [27] 
and  validated  somatic  variants  were  visually  curated 
using  IGV  [30]. 

Results 

Patient  selection  and  clinical  course 
We  selected  three  HGSC  patients  for  whom  matched 
ascites  tumour  specimens  were  available  from  the 
time  of  primary  surgery  and  two  subsequent  recur¬ 
rences.  As  fluid,  ascites  had  the  advantage  of  provid¬ 
ing  serial  samples  unrestricted  to  any  specific  tumour 
region.  As  shown  in  Figure  1,  patient  1  received  adju¬ 
vant  chemotherapy  with  five  cycles  of  carboplatin 
and  paclitaxel.  In  response  to  recurrence  of  disease 
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Years  from  diagnosis 

Figure  1.  Clinical  course  for  the  three  HGSC  patients  from  this 
study.  Black  squares  indicate  CA-1 25  levels  (kU/l)  over  time.  Arrows 
indicate  collection  of  ascites  tumour  specimens.  Black  bars  indicate 
chemotherapy.  Stars  indicate  time  of  death.  CATX,  carboplatin 
with  paclitaxel:  CA,  carboplatin;  TX,  taxol;  CAG,  carboplatin  with 
gemcitabine;  ETO,  etoposide;  GM,  gemcitabine;  DX,  pegylated 
liposomal  doxorubicin. 

8  months  after  completion  of  the  first  round  of  therapy, 
she  received  six  cycles  of  carboplatin  with  gemcitabine, 
which  slowed  tumour  progression  for  several  months. 
Upon  further  progression,  she  received  two  courses  of 
etoposide,  to  which  she  failed  to  respond.  Patient  2 
received  adjuvant  chemotherapy  with  six  cycles  of  car¬ 
boplatin.  In  response  to  progressive  disease,  she  was 
treated  with  seven  cycles  of  paclitaxel,  which  stabilized 
her  disease  for  about  2  months.  Upon  further  progres¬ 
sion,  she  received  one  cycle  of  gemcitabine,  to  which 
she  failed  to  respond.  Finally,  patient  3  received  adju¬ 
vant  chemotherapy  with  six  cycles  of  carboplatin  and 
paclitaxel,  which  led  to  complete  remission.  Upon  dis¬ 
ease  recurrence  at  9  months,  she  received  six  cycles 
of  carboplatin  with  gemcitabine,  which  led  to  a  sec¬ 
ond  complete  remission.  Upon  disease  recurrence,  she 
received  three  cycles  of  pegylated  liposomal  doxoru¬ 
bicin,  to  which  she  failed  to  respond. 

Identification  of  somatic  mutations  by  whole 
exome  sequencing 

Whole  exome  sequencing  was  performed  to  profile 
somatic  mutations  in  ascites  tumour  samples.  Single 
nucleotide  variants  and  insertions/deletions  (indels) 
were  validated  by  targeted  re-sequencing  at  a  depth 
of  >  50  000  reads  using  germline  DNA  for  comparison. 
We  identified  a  total  of  22-40  somatic  variants 
(somatic  point  mutations  plus  somatic  indels)  per 
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patient  (Tables  1  and  2).  Table  1  lists  all  somatic  point 
mutations,  their  allelic  frequencies,  and  their  predicted 
effects  on  protein  function  according  to  the  SIFT  pro¬ 
gram  [38].  Of  the  88  somatic  point  mutations  that  we 
found,  only  two  have  been  previously  reported,  one  in 
BRCA1  and  the  other  in  CSMD3  [39].  Of  the  remain¬ 
ing  86  somatic  point  mutations,  69  were  within  genes 
previously  found  to  be  mutated  in  ovarian  serous  carci¬ 
noma  according  to  the  Catalogue  of  Somatic  Mutations 
in  Cancer  (COSMIC)  database  [40].  By  SIFT  analysis, 
28%  of  somatic  point  mutations  were  predicted  at  high 
confidence  to  impair  protein  function.  In  addition  to 
somatic  point  mutations,  we  found  five  somatic  indels, 
which  involved  the  genes  TP  5  3  (patient  1),  MAC  FI 
(patient  2),  and  CAPN7,  DMD,  and  OR5A1  (patient  3) 
(Table  2).  The  only  gene  found  to  contain  a  somatic 
point  mutation  or  indel  in  all  three  patients  was  TP53. 

Clonal  structure  and  evolution  of  tumours 

Somatic  variants  were  subjected  to  high-depth  targeted 
re-sequencing  using  tumour  samples  from  all  three  time 
points.  This  allowed  quantification  of  allele  frequencies 
over  time,  even  for  mutations  that  fell  below  the  limit 
of  detection  in  the  original  exome  sequencing  libraries. 
We  found  that  the  vast  majority  of  mutations  that  were 
present  in  relapse  samples  were  present  in  primary 
samples  (Table  1).  Likewise,  the  majority  of  mutations 
that  were  present  in  primary  samples  persisted  in 
relapse  samples  (Table  1). 

We  used  a  hierarchical  clustering  model  to  infer 
tumour  clonality  based  on  concordant  changes  in 
mutant  allele  frequencies  (Figure  2).  Mutations  from 
patient  1  formed  three  distinct  clusters,  which  may 
reflect  three  distinct  clonal  populations.  Cluster  1  was 
markedly  more  abundant  than  the  other  two  clusters 
and  contained  a  heterozygous  TP53  mutation.  The 
relapse  tumours  from  patient  1  showed  a  doubling  in 
allelic  frequency  of  the  TP53  mutation,  suggesting  a 
loss  of  heterozygosity  at  the  TP 5 3  locus  (Table  1). 
The  relative  proportions  of  the  three  clusters  did  not 
change  significantly  across  the  three  time  points.  One 
exception  was  the  PLXB2  mutation  in  cluster  1,  which 
occurred  at  a  frequency  of  25%  in  the  primary  tumour 
sample  but  close  to  100%  in  relapse  samples.  This  sug¬ 
gests  that  the  PLXB2  mutation  occurred  after  the  other 
cluster  1  mutations  and  that  a  subclone  expressing  this 
mutation  increased  in  frequency  at  relapse  compared 
with  primary  tissue.  Patient  2  showed  four  major 
clusters  of  mutations.  As  with  patient  1,  the  cluster  con¬ 
taining  mutant  TP 5 3  was  dominant  at  all  time  points. 
Patient  3  showed  five  clusters  of  mutations  and  much 
more  complex  dynamics,  possibly  reflecting  the  dra¬ 
matic  changes  in  tumour  burden  seen  with  this  patient 
during  chemotherapy  (Figure  1).  Again,  the  cluster 
containing  mutant  TP 5 3  (cluster  1)  was  dominant  at 
all  time  points.  Cluster  3  increased  in  prevalence  at 
first  relapse  and  remained  prominent  at  second  relapse. 
Cluster  5  showed  the  opposite  pattern,  being  prevalent 
in  the  primary  tumour  but  absent  at  relapse. 
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Functional  annotation  of  mutant  genes 

Mutant  genes  from  each  cluster  were  functionally 
annotated  using  the  Gene  Ontology  (GO)  [35].  Despite 
the  lack  of  shared  mutations  among  patients,  the  muta¬ 
tions  fell  into  several  shared  functional  categories. 
The  predominant  cluster  from  all  three  patients  had 
cell  cycle  checkpoint  function.  In  addition,  statisti¬ 
cally  significant  functional  annotations  that  involved 
Golgi  vesicles  were  found  in  cluster  3  from  patient 
1  (FDR  =  9.3  x  10-14),  cluster  4  from  patient  2 
(FDR  =1.47  x  10-2),  and  cluster  2  from  patient  3 
(FDR  =1.38  x  1(T8). 

Analysis  of  shared  pathways  among  somatic 
mutations  and  germline  variants 

Germline  mutations  can  increase  the  risk  of  developing 
cancer  but  further  somatic  mutations  are  needed  to 
cause  the  actual  onset  of  the  disease.  In  HGSC, 
the  BRCA1,  BRCA2,  and  TP53  genes  frequently 
harbour  somatic  mutations,  and  germline  mutations 
can  contribute  to  a  familial  predisposition.  In  gen¬ 
eral,  tumourigenesis  and  disease  progression  require 
multiple,  parallel  pathways  that  can  be  altered  by 
either  germline  or  somatic  genetic  lesions  [41].  SNPs 
that  are  shared  among  these  patients  may  constitute 
a  common  predisposition  to  perturbed  pathways  that 
can  become  altered  further  by  somatic  mutations.  We 
asked  whether  in  these  patients,  somatically  mutated 
genes  were  related  to  genes  carrying  germline  coding 
variation.  We  identified  972  high-confidence  homozy¬ 
gous,  non-synonymous  SNPs  that  were  common  to 
all  three  patients.  These  SNPs  were  contained  within 
804  distinct  genes  (Supplementary  Table  1).  Of  the 
germline  SNPs  that  resulted  in  SIFT  scores,  3% 
were  predicted  at  high  confidence  to  impair  protein 
function.  Using  the  Cytoscape  GeneMania  plugin 
[32]  and  the  Reactome  FI  plugin  [37],  we  performed 
network  analysis  on  both  the  shared  homozygous 
germline  variants  and  the  somatic  variants  (Figure  3). 
The  network  revealed  connectivity  among  somatic 
and  germline  variants  that  defined  highly  connected 
modules.  The  top  two  interacting  nodes  that  were  rep¬ 
resented  as  network  hubs  were  TP53  and  fibronectin 
1 .  The  sets  of  mutated  genes  within  each  module  were 
then  assessed  for  functional  annotation  based  on  GO 
terms.  The  module  with  the  largest  number  of  nodes 
(57)  had  a  GO  annotation  of  ‘DNA  repair’  with  a 
g-value  of  2.8  x  10-8.  The  other  functionally  catego¬ 
rized  modules  were  ‘extracellular  matrix’  (40  nodes, 
g-value=1.0  x  10-20),  ‘G2/M  transition  of  mitotic 
cell  cycle’  (22  nodes,  q - value  =  3.3  x  10  7),  ‘trans- 
Golgi  network’  (13  nodes,  g-value  =  5.6  x  10-3), 
and  ‘Golgi  vesicle  targeting/membrane  budding’  (12 
nodes,  g-value=l.l  x  10-7).  Within  the  modules 
characterized  as  ‘DNA  repair’  and  ‘Golgi  vesicle 
targeting/membranebudding’,  there  was  enrichment  of 
somatic  mutations  relative  to  germline  variants  [11/57 
(19.3%)  and  2/7  (28.6%),  respectively]  when  compared 
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Table  1.  Validated  non-synonymous  somatic  mutations  ranked  by  allelic  frequencies  for  each  of  the  three  HGSC  patients 


Hg18  genomic 
position 

Gene  ID 

Type  of 
mutation 

Allele 
+  change 

Mutant  allele  frequency  by  tissue 

Amino  1st 

acid  Normal  Primary  relapse 

change  PBMC  ascites  ascites 

2nd 

relapse 

ascites 

SIFT  SIFT  dbSNP 

score  prediction*  ID 

COSMIC 

Patient  1 

Chrl  7:38497955 

BRCA1 

NS 

C>T 

S744N 

52-9 

65-1 

97-5 

96-7 

0.08 

T 

rs4986852:T 

Yes 

Chrl  7:35808624 

T0P2A 

NS 

T>C 

Y1163C 

0-1 

50-2 

97-1 

95-4 

0 

D 

Novel 

Yes 

Chr8:1 13310264 

CSMD3 

NS 

T>G 

N3621H 

49-7 

59-2 

93-0 

94-3 

0.56 

T 

rs  1 592624:G 

Yes 

Chr22:4907071 1 

PLXNB2 

NS 

C>T 

V144M 

0-3 

27-8 

92-8 

85-2 

0.01 

D 

Novel 

Yes 

Chrl  5:46367984 

SLC12A1 

NS 

G>A 

R951H 

0-1 

25-4 

56-3 

64-1 

N/A 

NS 

Novel 

Yes 

Chr2:136589859 

CXCR4 

US 

T>A 

N41Y 

0-1 

15-2 

47-9 

56-2 

N/A 

NS 

Novel 

No 

Chr10:92972714 

PCGF5 

NS 

C>A 

H36N 

0-1 

23-8 

40-4 

53-2 

0 

D 

Novel 

No 

Chr19:61061829 

NLRP4 

NS 

C>T 

R420W 

0-2 

25-9 

57-7 

47-2 

0 

D 

Novel 

Yes 

Chr8 : 1 41 637781 

EIF2C2 

NS 

G>T 

P288H 

0-1 

29-7 

37-1 

44-7 

0 

D 

Novel 

No 

Chrl  5:31 659551 

RYR3 

NS 

C>T 

Q451* 

0-2 

18-4 

25-1 

35-4 

N/A 

N/A 

Novel 

Yes 

Chr19:6616095 

TNFSF14 

NS 

G>A 

R189W 

0-1 

17-7 

30-1 

30-5 

0 

DLC 

Novel 

Yes 

Chrl  4:1 9834390 

TTC5 

NS 

G>A 

T232M 

0-4 

12-9 

40-5 

29-7 

0.14 

T 

Novel 

Yes 

Chr8:691 19340 

PREX2 

NS 

A>C 

T302P 

0-0 

0-0 

19-1 

28-6 

0 

DLC 

Novel 

Yes 

Chrl  4:41 430594 

LRFN5 

NS 

G>T 

A593S 

0-1 

15-8 

16-5 

28-2 

0.23 

T 

Novel 

Yes 

Chr20:4798604 

SLC23A2 

NS 

C>T 

A400T 

0-1 

22-3 

41-6 

27-7 

0.24 

T 

Novel 

Yes 

Chrl  7:5374631 

NLRP1 

NS 

C>G 

E1142D 

0-2 

15-7 

17-6 

26-7 

0.3 

T 

Novel 

Yes 

Chr3:123608886 

F AMI  62 A 

NS 

C>T 

T111M 

0-1 

9-1 

20-3 

26-2 

0 

D 

Novel 

Yes 

Chr3:1 741 57167 

SPATA16 

NS 

A>G 

Y359H 

0-1 

4-8 

28-8 

23-1 

0 

DLC 

Novel 

Yes 

Chr9:140135932 

CACNA1B 

NS 

G>A 

G2165S 

0-1 

20-4 

28-1 

22-6 

0.76 

T 

Novel 

Yes 

Chrl  4:94023488 

SERPINA12 

NS 

C>T 

V384I 

0-0 

9-7 

12-4 

21-8 

0.95 

T 

Novel 

Yes 

Chr3:1 40576092 

C0PB2 

NS 

C>T 

G227E 

0-0 

12-3 

27-3 

17-9 

0 

D 

Novel 

Yes 

Patient  2 

Chr8:9481 5372 

RBM12B 

US 

G>A 

P815S 

50-0 

68-0 

97-9 

99-9 

N/A 

NS 

Novel 

Yes 

Chrl  7:751 8901 

TP53 

PSS 

C>T 

N/A 

0-4 

99-7 

99-1 

99-8 

N/A 

N/A 

Novel 

Yes 

Chr19:3498535 

C19orf28 

NS 

A>C 

L283R 

47-4 

99-8 

99-8 

99-8 

0 

DLC 

Novel 

No 

ChrX:1 2551 3552 

DCAF12L1 

NS 

C>G 

V241L 

0-2 

99-1 

91-7 

98-9 

0.15 

T 

Novel 

Yes 

Chr21 :4681 1715 

DIP2A 

US 

G>T 

V1447L 

0-7 

66-4 

66-4 

78-1 

N/A 

NS 

Novel 

Yes 

Chr10:120791840 

EIF3A 

NS 

C>T 

R1061H 

0-4 

39-5 

62-6 

62-2 

0.02 

DLC 

Novel 

No 

Chr19:47483139 

CIC 

NS 

G>C 

G120A 

0-5 

51-3 

59-5 

55-4 

0.65 

T 

Novel 

Yes 

Chrl  1:14810942 

PDE3B 

NS 

T>A 

Y731  * 

0-6 

52-4 

45 

53-8 

N/A 

N/A 

Novel 

No 

Chr2:235069395 

ARL4C 

NS 

C>G 

R 1 92  P 

0-2 

38-2 

36-0 

53-3 

0 

DLC 

Novel 

Yes 

Chr9:1 30385899 

SPTAN1 

NS 

C>T 

R675C 

0-2 

41-7 

56-7 

51-6 

0.01 

DLC 

Novel 

Yes 

Chr4:187779522 

FAT1 

US 

G>A 

Q1741* 

0-3 

29-8 

47 

51-1 

N/A 

NS 

Novel 

Yes 

Chr9:1 30061 400 

G0LGA2 

NS 

T>C 

Q616R 

0-5 

46-7 

42-6 

50-5 

0.4 

T 

Novel 

Yes 

Chr5:71775527 

ZNF366 

NS 

G>T 

L683I 

0-9 

45-1 

54-2 

49-3 

0.32 

T 

Novel 

No 

Chr9:91 198306 

SEMA4D 

NS 

C>T 

VI 671 

0-3 

64-6 

57-6 

47-7 

0.23 

T 

Novel 

Yes 

Chrl  1:51 268538 

0R4A5 

NS 

A>G 

V145A 

0-1 

30-7 

41-1 

40-2 

0.33 

T 

Novel 

Yes 

Chr17:33738996 

G  PR  179 

US 

G>C 

R1328G 

0-5 

21-4 

48-1 

37-5 

N/A 

NS 

Novel 

Yes 

Chr7:1 47545382 

CNTNAP2 

NS 

G>A 

R1027K 

0-1 

15-9 

34-5 

37-3 

1 

T 

Novel 

Yes 

ChrX:22005541 

PHEX 

NS 

C>G 

P155A 

0-0 

0-0 

0-1 

35-5 

0 

D 

Novel 

Yes 

Chrl  :31923029 

C0L16A1 

NS 

G>A 

R677C 

0-3 

47-9 

35-3 

34-2 

0.01 

D 

Novel 

Yes 

Chr10:104818466 

CNNM2 

NS 

C>A 

P723T 

0-4 

27-2 

31-5 

34-1 

0.52 

T 

Novel 

Yes 

Chr2:21085858 

APOB 

NS 

C>A 

A2463S 

0-6 

56-9 

31-1 

27-6 

0.09 

T 

Novel 

Yes 

Chrl  6:20556245 

ACSM1 

NS 

C>T 

Mil 

0-1 

0-1 

0-1 

21-5 

0.02 

DLC 

Novel 

No 

Chrl  :32034882 

SP0CD1 

NS 

G>A 

Q21 6* 

0-1 

0-1 

0-0 

19-7 

N/A 

N/A 

Novel 

Yes 

Chr22:28075220 

AP1B1 

NS 

T>C 

D475G 

0-0 

0-0 

0-0 

19-3 

0.02 

D 

Novel 

Yes 

Chrl  :1 48124765 

HIST2H2BE 

NS 

T>C 

K17R 

0-3 

22-3 

21-3 

18-3 

0.05 

DLC 

Novel 

Yes 

Chr20:32934365 

ACSS2 

NS 

A>T 

M96L 

0-4 

16-6 

24-6 

11-0 

0 

D 

Novel 

No 

Chr3:13345395 

NUP210 

NS 

G>C 

Q1388E 

0-1 

39-8 

39-2 

9-3 

N/A 

NS 

Novel 

Yes 

Chr10:95921058 

PLCE1 

NS 

C>T 

Q542* 

0-0 

0-0 

0-1 

8-4 

N/A 

N/A 

Novel 

Yes 

Chrl  6:2881 6828 

ATP2A1 

NS 

C>A 

R476S 

0-1 

4-5 

19-9 

1-4 

0.4 

T 

Novel 

Yes 

Chr8:88954771 

DCAF4L2 

NS 

A>G 

I182T 

0-3 

24-1 

0-1 

0-0 

0.01 

DLC 

Novel 

Yes 

Patient  3 

Chr6:1 46761 919 

GRM1 

NS 

G>T 

R684L 

0-1 

83-9 

65-3 

96-7 

0 

D 

Novel 

Yes 

Chrl  8:75996662 

ADNP2 

NS 

G>C 

R792T 

0-0 

86-2 

43-6 

96-1 

0.48 

T 

Novel 

Yes 

Chr14:73894317 

Cl4orf115 

NS 

T>A 

C360S 

0-0 

93-9 

65-6 

94-2 

0.95 

T 

Novel 

No 

Chrl  7:751 901 5 

TP53 

PSS 

C>T 

N/A 

0-0 

92-5 

74-6 

93-5 

N/A 

N/A 

Novel 

Yes 

ChrX:47200359 

ZNF41 

NS 

C>T 

G67E 

0-1 

56-4 

54-1 

69-6 

0 

D 

Novel 

Yes 

ChrX:1 3525461 9 

GPR112 

NS 

A>G 

E300G 

0-1 

46-9 

54-6 

67-0 

0 

DLC 

Novel 

Yes 

Chrl  6:1 1839419 

RSL1D1 

US 

C>G 

S400T 

0-0 

58-0 

46-5 

64-8 

N/A 

NS 

Novel 

No 

Chr4:1 201 70760 

5YNP02 

NS 

A>C 

N461T 

0-0 

45-7 

52-8 

62-9 

0.31 

T 

Novel 

Yes 

Chrl  6:82722355 

HSDL1 

NS 

G>C 

L25V 

0-0 

3-5 

55-0 

60-2 

0 

D 

Novel 

No 

Chrl  7:1 6553073 

CCDC144A 

NS 

G>C 

R326T 

0-0 

0-0 

35-2 

60-0 

0 

DLC 

Novel 

No 

Chrl  1:88794758 

N0X4 

NS 

A>T 

FI  971 

0-0 

20-5 

45-0 

56-5 

0 

D 

Novel 

Yes 

Copyright  ©  2012  Pathological  Society  of  Great  Britain  and  Ireland. 
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Table  1.  (Continued) 


Hg  1 8  genomic 
position 

Gene  ID 

Type  of 
mutation+ 

Allele 

change 

Mutant  allele  frequency  by  tissue 

Amino  1st 

acid  Normal  Primary  relapse 

change  PBMC  ascites  ascites 

2nd 

relapse 

ascites 

SIFT 

score 

SIFT 

prediction* 

dbSNP 

ID 

COSMIC 

Chr3:1 30452227 

COPG 

NS 

A>G 

N17S 

0.0 

50-1 

38-3 

55-0 

0.32 

T 

Novel 

Yes 

Chrl  2 :48 141 063 

SPATS2 

NS 

A>G 

M1V 

0.0 

51-3 

51-2 

52-9 

0.29 

T 

Novel 

No 

Chr9:108728098 

ZNF462 

NS 

G>A 

R695K 

0.0 

44-1 

37-6 

52-2 

0 

DLC 

Novel 

Yes 

Chrl  2 : 1 038043 1 0 

SLC41A2 

NS 

C>G 

K223N 

0.1 

22-8 

23-1 

45-8 

0.06 

T 

Novel 

Yes 

Chr9:73549906 

TMEM2 

NS 

G>T 

S294R 

0.0 

31-3 

24-1 

40-1 

0 

D 

Novel 

Yes 

Chr19:45234235 

ZNF780B 

US 

T>C 

Y168C 

0.0 

2-1 

26-1 

40-0 

N/A 

NS 

Novel 

Yes 

Chrl  1:551 27965 

0R4C1 1 

NS 

G>T 

S154Y 

0.0 

0-5 

26-2 

35-5 

0 

D 

Novel 

Yes 

Chr10:1 25629803 

CPXM2 

NS 

T>C 

N106S 

0.1 

47-1 

28-8 

33-2 

0.14 

T 

Novel 

Yes 

Chr2:197302318 

CCDC150 

NS 

G>C 

A392P 

0.0 

27-2 

31-3 

32-9 

0.16 

T 

Novel 

No 

Chrl  2:69256442 

PTPRB 

US 

G>C 

S725R 

0.1 

2 

26-2 

31-9 

N/A 

NS 

Novel 

No 

Chrl  0:1 0267381 8 

FAM178A 

NS 

G>T 

R357I 

0.0 

23-6 

32-9 

30-0 

0.01 

DLC 

Novel 

No 

Chrl  :152190704 

CRTC2 

NS 

C>A 

A34S 

0.0 

34-3 

23-9 

29-7 

0 

DLC 

Novel 

Yes 

Chrl  :165933526 

RCSD1 

NS 

G>C 

E347D 

0.0 

24-8 

19 

23-4 

0 

DLC 

Novel 

Yes 

Chrl  2:64990591 

HELB 

NS 

C>G 

P539R 

0.0 

0-1 

25-9 

22-0 

0 

D 

Novel 

No 

Chr3:9966618 

PRRT3 

NS 

A>T 

F61Y 

0.1 

39-9 

19 

20-9 

0 

DLC 

Novel 

No 

Chr6:571 20371 

ZNF451 

NS 

C>G 

S510* 

0.0 

0-6 

19-6 

20-8 

N/A 

N/A 

Novel 

Yes 

Chrl  :6627704 

DNAJC1 1 

NS 

C>G 

V322L 

0.0 

2-6 

24 

20-1 

1 

T 

Novel 

Yes 

Chrl  :1 84381 645 

HMCN1 

NS 

C>G 

Q4859E 

0.0 

19-8 

18-1 

16-8 

0 

DLC 

Novel 

Yes 

Chrl  5:1 9884381 

0R4N4 

NS 

G>A 

R182Q 

0.1 

1-8 

12-9 

16 

0.28 

T 

Novel 

Yes 

Chrl  2:877064 

WNK1 

NS 

C>G 

S1889C 

0.0 

9-4 

0-0 

0-0 

0 

DLC 

Novel 

Yes 

Chr5:1 24010504 

ZNF608 

NS 

A>T 

S1158T 

0.1 

18-1 

0-1 

0-0 

0.16 

T 

Novel 

Yes 

Chr12:13099988 

KIAA1467 

NS 

C>G 

Q92E 

0.0 

27-2 

0-0 

0-0 

1 

T 

Novel 

Yes 

Chrl  1:63419288 

MARK2 

NS 

G>A 

E13K 

0.0 

28-9 

0-0 

0-0 

0.08 

T 

Novel 

Yes 

Chrl  4:72072653 

RGS6 

NS 

A>C 

I429L 

0.0 

44-4 

0-0 

0-0 

0 

D 

Novel 

No 

Chrl  7:57022774 

NACA2 

NS 

C>G 

V184L 

0.1 

26-6 

0-1 

0-0 

0 

D 

Novel 

Yes 

Chr12:38234148 

ABCD2 

NS 

C>G 

E686Q 

0.0 

16-9 

0-0 

0-0 

0.44 

T 

Novel 

Yes 

*Stop  codon. 

+  NS  =  non-synonymous  substitution;  US  =  unknown  substitution;  PSS=  possible  splice  site. 

*D  =  damaging;  DLC  =  damaging  (low  confidence);  T  =  tolerated;  NS  =  not  scored;  N/A=  not  available. 


Table  2.  Validated  indels  identified  in  the  three  HGSC  patients 


Hg  1 8  genomic 
position 

Gene  ID 

Type  of 
mutation 

Allele  change 

COSMIC 

mutation 

dbSNP  ID 

Origin 

of  indel  by  patient 

Patient  1 

Patient  2 

Patient  3 

Chr14:22618623 

ACINI 

Insertion 

AG AACGTG AACGT G A> 

Yes 

rs  14803 158 

Germline 

Germline 

AG AACGT G AACGT  - 

G AACGTG A 

Chr7:1 50344835 

ATG9B 

Insertion 

GC>GCC 

Yes 

novel 

Germline 

Germline 

Chrl  1:76429189 

B3GNT6 

Deletion 

CTT>CT 

No 

novel 

Germline 

Germline 

Germline 

Chrl  9:40950777 

Cl9orf55 

Deletion 

AGGG>AGG 

No 

rs66706567 

Germline 

Germline 

Germline 

Chr3:1 5263894 

CAPN7 

Deletion 

T  C  AC  A  AAA  AT  A  ACCCC  A  >  T  CA 

Yes 

novel 

Somatic 

ChrX:3 1435324 

DMD 

Deletion 

GTT  G  AG  AG  ACTTTT  >  GTTT 

Yes 

novel 

Somatic 

Chr2:1 76696536 

H0XD9 

Insertion 

CG  C AG  C AG  C AG  C A  > 

No 

rs66706567 

Germline 

Germline 

CGCAGCAGCAGCAGCA 

Chr14:93651883 

IFI27 

Insertion 

TGCC>TGGCCATGGCGGC 

No 

novel 

Germline 

Germline 

Germline 

Chrl  :39553852 

MACF1 

Deletion 

AGT CT CCAT CT AGTT CAAGT GTCC >  Yes 

novel 

Somatic 

AGTCTAGTGTCC 

Chrl  1:58967477 

0R5A1 

Deletion 

ACTT  CTT  CT  >  ACTT  CT 

Yes 

novel 

Somatic 

Chrl  7:3541 025 

P2RX5/TAX1BP3 

Deletion 

TGGGGG>TGGGG 

Yes 

rsl  50260459 

Germline 

Germline 

Chr2:1 13710767 

Pax8 

Insertion 

GCC>GCC 

No 

rs5826760 

Germline 

Germline 

Germline 

Chrl  9 :229 1  154 

SPPL2B 

Insertion 

GCC>GCCC 

No 

novel 

Germline 

Germline 

Germline 

Chr7:149107599 

SSPO 

Insertion 

GCCCC>GCCCCC 

No 

novel 

Germline 

Germline 

Germline 

Chr20:421 27972 

T0X2 

Deletion 

CCCGCCGCCGCC> 

Yes 

novel 

Germline 

CCCGCCGCCGCCGCC 

Chrl  7 :75201 14 

TP53 

Deletion 

GGGAAGG>GGG 

Yes 

novel 

Somatic 

Chr4:75995 

ZNF595 

Insertion 

CA>CAA 

No 

rsl  39953554 

Germline 

Germline 

Germline 
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BRCA1 
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TP53 


SLC12A1 
NLRP4 
CXCR4 
PCGF5 
EIF2C2 
RYR3 
TNFSF14 
CACNA1 B 
TTC5 
SLC23A2 
PREX2 
SPATA16 
LRFN5 
NLRP1 
FAM162A 
SERPINA12  | 
C0PB2 


RBM12B 

DIP2A 

TP53 

C19orf28 

DCAF12L1 

EIF3A 

CIC 

SPTAN1 

SEMA4D 

PDE3B 

G0LGA2 

ARL4C 

FAT1 

OR4A5 

GPR179 

CNTNAP2 

CNNM2 

MAC  FI 

C0L16A1 

APOB 

NUP210 

PHEX 

ACSM1 

SPOCD1 

AP1B1 

PLCE1 

HIST2H2BE 

ACSS2 

ATP2A1 

DCAF4L2 


Top  GO  annotation 

DNA  repair  (fdr  i.isx  io'5) 
cell  cycle  checkpoint 

(FDR  2.58  x  1 0'2) 

axon  guidance 

(FDR  1.15  x  10'2) 

gene  silencing  by  RNA 

(FDR  8.05  x  1 0'11) 


CO  PI  coating  of 
Golgi  vesicle 

(FDR  9.3  x  10'14) 

divalent  metal  ion 
transport 

(FDR  1.37  x  1 0'9) 


signal  transduction 
involved  in  mitotic  cell 
cycle  checkpoint 

(FDR  1.68  x  1 0'1) 


eukaryotic  translation 
initiation  factor  3 
complex 

(FDR  1.55  x  1 0'7) 


nuclear-transcribed 
mRNA  catabolic  process, 
exonucleolytic 

(FDR  4.77  x  1 0'9) 


post-Golgi  vesicle- 
mediated  transport 

(FDR  1.47  x  1 0'2) 

canonical  Wnt  receptor 
signalling  pathway 

(FDR  1.47  x  1 0'2) 


signal  transduction  involved  in 
mitotic  cell  cycle  checkpoint 

(FDR  1.68  x  1 0'1) 

COPI  coating  of  Golgi  vesicle 

(FDR  1.38  x  1 0'8) 

dystrophin-associated 
glycoprotein  complex 

(FDR  3.96x1  O'6) 


replisome 

(FDR  7.88  x  1 0'1) 


protein  serine/threonine 
kinase  activity 

(FDR  1.78  x  1 0'1) 


microbody  membrane 

(FDR  2.07  x  1 0'3) 

peroxisome  membrane 

(FDR  2.07  x  1 0'3) 


Per  cent  Frequency  Scale 


Patient  1  Patient  2  Patient  3 


34  31  40 


94  75  97 


Figure  2.  Mutant  allele  frequency  trend  and  inference  of  tumour  clonality.  Mutant  allele  frequencies  (0-100%)  are  represented  by  the 
green-red  heat  map,  according  to  the  indicated  scale.  Clonality  was  inferred  using  a  hierarchical  Ward  clustering  algorithm  based  on 
temporal  changes  in  non-synonymous  mutant  allele  frequencies.  Gene  Ontology  (GO)  annotation  of  function  was  performed  on  the  genes 
specified  by  each  cluster  using  the  Cytoscape  Genemania  plugin.  (A)  Patient  1 ;  (B)  patient  2;  (C)  patient  3. 
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Figure  3.  Functional  interactions  network  analysis  of  germline  and  somatic  non-synonymous  coding  variants.  Direct  interactions  are 
illustrated  by  connecting  edges  and  node  size  is  relative  to  the  number  of  interactions.  Genes  that  contained  shared,  homozygous,  germline 
variants  (circular  nodes)  in  all  three  patients  and  somatic  mutations  (diamond  nodes)  were  analysed  by  the  network  analysis  tool  Cytoscape 
Genemania  plugin.  Genes  predicted  by  GeneMania  are  shown  as  hexagonally  shaped  nodes.  Data  sources  used  for  the  predictions  included 
'pathway'  data,  which  describe  two  genes  that  share  a  reaction,  or  'physical'  data,  which  describe  protein -protein  interactions.  Node 
colours  denote  genes  that  belong  to  network  modules  with  annotated  GO  terms. 


with  the  total  number  of  somatic  versus  germline  vari¬ 
ants,  93/804  (11.6%).  The  modules  for  ‘extracellular 
matrix’  and  ‘G2/M  transition  of  mitotic  cell  cycle’  also 
contained  both  somatic  mutations  and  germline  vari¬ 
ants.  These  results  suggest  possible  interplay  between 
SNPs  that  may  impart  a  hereditary  predisposition  to 
HGSC  and  subsequently  acquired  somatic  mutations. 


Discussion 

Epithelial  ovarian  cancer  is  highly  responsive  to  front¬ 
line  treatment,  but  relapse  with  chemotherapy-resistant 
disease  occurs  in  the  majority  of  cases.  To  assess 
disease  progression,  we  used  CD45—  sorted  ascites 
cells  as  a  surrogate  for  solid  tumour  biopsies  since 
recurrent  samples  are  not  surgically  resected  and 
ascites  is  the  only  available  source  of  tumour  in  most 
cases.  We  sequenced  the  exomes  of  primary,  first, 
and  second  relapse  tumours  and  discovered  muta¬ 
tions  affecting  multiple  cellular  processes  and  struc¬ 
tures,  including  cell  cycle  control,  Golgi  vesicles,  DNA 
repair,  and  extracellular  matrix.  Remarkably,  89%  of 
mutations  found  in  relapse  tumours  were  present  in 
matched  primary  tumours,  indicating  that  recurrent 
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HGSC  arises  from  the  selective  persistence  or  out¬ 
growth  of  pre-existing  tumour  clones,  with  relatively 
little  accumulation  of  new  mutations.  Our  results  raise 
the  encouraging  possibility  that  improved  targeting  of 
mutated  pathways  in  primary  tumours  could  potentially 
prevent  the  development  of  recurrent  disease. 

Our  finding  that  the  majority  of  mutations  in  recur¬ 
rent  disease  were  present  in  the  primary  tumour  sam¬ 
ple  is  consistent  with  a  prior  study  of  HGSC  that 
used  FISH,  SNP  arrays,  and  array  CGH  [24].  Given 
the  genotoxic  effects  of  platinum-based  chemother¬ 
apy,  it  is  surprising  that  so  few  new  mutations  were 
found  in  relapse  tumours.  Rather,  our  data  indicate 
that  chemotherapy  failed  to  destroy  the  major  clones 
present  in  the  primary  tumour,  as  indicated  by  the 
fact  that  fewer  than  10%  of  mutations  were  lost  after 
treatment.  An  exception  was  patient  3,  who  responded 
dramatically  to  her  first  round  of  treatment  based 
on  decreased  CA-125  levels  (Figure  1)  and  complete 
disappearance  of  a  mutational  cluster  (cluster  5  in 
Figure  2).  Patient  3  also  showed  the  greatest  num¬ 
ber  of  new  mutations  in  relapse  samples.  Relapse- 
specific  mutations  observed  in  patients  1  and  2  included 
PREX2 ,  which  has  recently  been  reported  in  metastatic 
melanomas  [42],  and  APlbl ,  which  is  associated  with 
vesicle  formation  in  the  trans-Golgi  network  [43].  It 
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is  unknown  whether  these  mutations  were  due  to  the 
DNA-damaging  effects  of  chemotherapy  or  the  expan¬ 
sion  of  a  minor  tumour  clone  from  the  primary  tumour. 

In  all  three  patients,  GO  annotated  somatic  muta¬ 
tion  clusters  were  seen  that  involved  either  the  Golgi 
vesicle  or  Golgi  transport  (Figure  2).  Similarly,  net¬ 
work  analysis  of  somatic  mutations  and  germline  SNPs 
revealed  two  modules  that  were  annotated  for  the  trans- 
Golgi  network  and  Golgi  vesicle  targeting  (Figure  3). 
The  Golgi  apparatus  has  been  previously  implicated 
in  cancer  processes  such  as  drug  efflux,  apoptosis, 
and  proliferation.  Cisplatin  is  sequestered  by  the  Golgi 
in  ovarian  carcinoma  cells,  which  may  be  necessary 
for  platinum  efflux  [44].  The  Golgi  has  also  been 
implicated  in  apoptosis  by  such  means  as  compartmen¬ 
talizing  caspase  2,  triggering  mitochondrial  membrane 
permeability,  releasing  death  receptors  to  the  cell  sur¬ 
face,  and  sensing  misfolded  proteins  [45].  Additionally, 
anti-cancer  drugs  such  as  2-phenyl  benzimidazole  (2- 
PB)  disrupt  the  Golgi  apparatus,  which  leads  to  reduced 
cell  proliferation  and  tumour  growth  [46]. 

Germline  SNPs  and  somatic  mutations  were  also 
found  in  the  extracellular  matrix  (ECM)  pathway, 
including  collagens,  laminins,  integrins,  syndecan  1, 
and  fibronectin  1.  These  proteins  are  involved  in  cell 
adhesion,  cell  signalling,  and  cytoskeletal  organization, 
and  alterations  to  ECM  pathways  have  been  found 
during  the  evolution  of  other  cancers  [20,47]  but  not 
previously  in  HGSC. 

Germline  and  somatic  variants  were  found  in  path¬ 
ways  such  as  p53,  cell  cycle  progression,  and  DNA 
repair  (Figure  3),  which  are  known  HGSC  suscepti¬ 
bility  pathways  [7].  TP 5 3  was  the  only  gene  found 
to  harbour  a  somatic  mutation  in  all  three  patients, 
while  BRCA1  was  somatically  mutated  in  patient  1 
and  all  three  patients  shared  the  same  germline  SNPs  in 
BRCA2  and  ATM  (Table  1,  Supplementary  Table  1  and 
Supplementary  Table  2).  An  apparent  chemotherapy- 
susceptible  clone,  cluster  5  from  patient  3,  contained 
a  highly  abundant  RGS6  mutation,  which  is  a  gene 
that  has  been  shown  to  be  associated  with  cancer 
risk  [48].  More  recently,  RGS6  has  been  reported  to 
activate  the  intrinsic  pathway  of  apoptosis  in  a  p53- 
independent  manner  and  affected  doxorubicin  suscep¬ 
tibility  in  breast  cancer  cells  [49]. 

Cancer  pathogenesis  can  occur  through  the  loss  of 
both  alleles  in  a  tumour  suppressor  gene  in  the  clas¬ 
sical  ‘two-hit’  hypothesis  of  tumour  suppression  or 
through  a  more  subtle,  titrated  loss  of  tumour  sup¬ 
pressor  gene  function  as  described  by  a  continuum 
model  of  tumour  suppression  [50].  Consistent  with 
the  continuum  model,  tumourigenesis  can  be  initi¬ 
ated  by  p53  haploinsufficiency,  as  seen  in  patients 
with  Li-Fraumeni  syndrome  and,  it  appears,  patient 
1  (Table  1).  The  continuum  model  of  disease  pro¬ 
gression  may  extend  beyond  tumour  suppressor  genes, 
since  we  observed  several  examples  of  pathways  that 
were  altered  by  germline  variants  and  became  further 
perturbed  by  somatic  mutations  (Figure  3). 
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We  have  shown  here  that  HGSC  tumours  are  het¬ 
erogeneous,  with  multiple,  genetically  distinct  clones 
present  prior  to  treatment.  Most  of  the  mutations  in 
primary  tumours  persist  despite  treatment,  suggesting 
that  most  clones  are  able  to  evade  current  chemother¬ 
apy  treatments.  Furthermore,  there  is  little  accumu¬ 
lation  of  new  mutations  in  recurrent  tumours,  even 
after  two  rounds  of  chemotherapy.  Encouragingly,  we 
found  novel  mutations  in  the  Golgi  and  ECM  path¬ 
ways,  which  may  offer  new  therapeutic  targets  for 
HGSC.  Moreover,  the  relatively  stable  nature  of  the 
HGSC  genome  over  time  provides  hope  that  improved 
therapeutic  targeting  of  the  genetic  lesions  present  in 
primary  tumours  will  ultimately  be  effective  for  pre¬ 
venting  tumour  recurrence. 
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Abstract 

Purpose:  Tumor-infiltrating  lymphocytes  (TIL),  in  particular  CD8+  T  cells  and  CD20+  B  cells,  are 
strongly  associated  with  survival  in  ovarian  cancer  and  other  carcinomas.  Although  CD8+  TIL  can  mediate 
direct  cytolytic  activity  against  tumors,  the  role  of  CD20+  TIL  is  poorly  understood.  Here,  we  investigate  the 
possible  contributions  of  CD20+  TIL  to  humoral  and  cellular  tumor  immunity. 

Experimental  Design:  Tumor  and  serum  specimens  were  obtained  from  patients  with  high-grade  serous 
ovarian  cancer.  CD8+  and  CD20+  TIL  were  analyzed  by  immunohistochemistry  and  flow  cytometry. 
Immunoglobulin  molecules  were  evaluated  by  DNA  sequencing.  Serum  autoantibody  responses  to  the 
tumor  antigens  p53  and  NY-ESO-1  were  measured  by  ELISA. 

Results:  The  vast  majority  of  CD20+  TIL  were  antigen  experienced,  as  evidenced  by  class-switching, 
somatic  hypermutation,  and  oligoclonality,  yet  they  failed  to  express  the  canonical  memory  marker  CD27. 
CD20+  TIL  showed  no  correlation  with  serum  autoantibodies  to  p53  or  NY-ESO- 1 .  Instead,  they  colocalized 
with  activated  CD8+  TIL  and  expressed  markers  of  antigen  presentation,  including  MHC  class  I,  MHC  class 
II,  CD40,  CD80,  and  CD86.  The  presence  of  both  CD20+  and  CD8+  TIL  correlated  with  increased  patient 
survival  compared  with  CD8+  TIL  alone. 

Conclusions:  In  high-grade  serous  ovarian  tumors,  CD20+  TIL  have  an  antigen- experienced  but  atypical 
CD27-  memory  B-cell  phenotype.  They  are  uncoupled  from  serum  autoantibodies,  express  markers  of 
antigen-presenting  cells,  and  colocalize  with  CD8+  T  cells.  We  propose  that  the  association  between  CD20+ 
TIL  and  patient  survival  may  reflect  a  supportive  role  in  cytolytic  immune  responses.  Clin  Cancer  Res;  1-12. 
©2012  AACR. 


Introduction 

High-grade  serous  ovarian  cancer  (hereafter  abbreviated 
HGSC)  is  the  most  common  and  lethal  subtype  of  ovarian 
cancer,  with  a  5 -year  survival  rate  of  only  30%  ( 1 ) .  Although 
most  patients  are  highly  responsive  to  primary  surgery  and 
chemotherapy,  the  majority  experience  recurrence  within 
1  to  3  years  and  ultimately  succumb  to  their  disease  (2). 
Despite  these  unfortunate  statistics,  a  subset  of  patients 
experience  prolonged  disease-free  survival,  favorable  prog- 
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nostic  indicators  for  HGSC  include  early  stage  and  optimal 
surgical  debulking  (2).  furthermore,  studies  over  the  past  2 
decades  have  revealed  strong  links  between  host  immunity 
and  survival.  In  particular,  multiple  studies  have  shown  that 
the  presence  of  CD3+  and  CD8+  tumor-infiltrating  lym¬ 
phocytes  (TIL)  is  associated  with  markedly  prolonged  sur¬ 
vival  (3-8). 

In  addition  to  CD8+  T  cells,  we  have  shown  that  tumor- 
infiltrating  CD20+  B  cells  (CD20+  TIL)  are  strongly  asso¬ 
ciated  with  patient  survival  in  HGSC  (6).  Tumor-infiltrating 
B  cells  are  also  correlated  with  favorable  outcomes  in  breast, 
cervical,  and  non-small  cell  lung  cancer  (9-11).  In  breast 
cancer  and  germ  cell  tumors,  tumor-infiltrating  B  cells  have 
been  shown  to  consist  of  activated,  antigen-experienced, 
oligoclonal  subpopulations  (12-18).  Beyond  this  informa¬ 
tion,  the  mechanistic  role  of  CD20+  TIL  in  tumor  immunity 
remains  poorly  defined.  A  priori,  tumor-infiltrating  B  cells 
could  mediate  their  effects  through  the  production  of 
tumor-specific  antibodies.  Indeed,  more  than  40%  of  HGSC 
patients  show  autoantibody  responses  to  tumor  antigens 
such  as  p53  and  NY-ESO-1  (19,  20).  In  theory,  tumor- 
specific  autoantibodies  could  mediate  antitumor  responses 
by  direct  inhibition  of  target  proteins,  activation  of 
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Translational  Relevance 

High-grade  serous  ovarian  cancer  has  a  5 -year  survival 
rate  of  only  30%,  so  there  is  an  urgent  need  for  improved 
treatments.  Recent  studies  have  provided  unequivocal 
evidence  that  patients  who  mount  strong  immune 
responses  against  their  tumor  experience  markedly  pro¬ 
longed  survival.  Although  tumor-infiltrating  T  cells  are 
clearly  important  for  tumor  immunity,  we  have  recently 
shown  an  equally  prominent  role  for  B  cells.  Here,  we 
investigate  the  underlying  mechanisms.  Unexpectedly, 
we  found  that  tumor-infiltrating  B  cells  were  uncoupled 
from  antibody  responses  to  common  tumor  antigens. 
Instead,  they  colocalized  with  T  cells  and  showed  char¬ 
acteristics  of  antigen-presenting  cells.  The  presence  of 
both  B  cells  and  T  cells  in  tumors  was  associated  with 
better  prognosis  than  T  cells  alone.  Thus,  tumor-infil¬ 
trating  B  cells  seem  to  support  T-cell  responses  to  cancer. 
Immunotherapeutic  strategies  that  engage  both  lympho¬ 
cyte  subsets  may  have  more  potent  and  sustained  anti¬ 
tumor  effects. 


complement,  opsonization  of  tumor  antigens,  or  triggering 
of  antibody-dependent  cellular  cytotoxicity.  In  addition  to 
secreting  antibodies,  B  cells  could  enhance  tumor  immunity 
through  cellular  mechanisms.  For  example,  B  cells  can  serve 
as  antigen-presenting  cells  (APC)  to  T  cells  (reviewed  in 
ref.  21).  In  this  capacity,  B  cells  have  an  advantage  over 
dendritic  cells  in  that  they  can  concentrate  their  cognate 
antigen  through  surface  Ig  molecules,  thereby  allowing 
presentation  of  even  low  abundance  antigens  (21). 

Additional  insights  into  the  possible  role(s)  of  CD20+  TIL 
in  tumor  immunity  can  be  gained  from  studies  of  B  cell 
infiltrates  in  autoimmunity  and  solid  organ  transplantation 
(21-29).  In  renal  allografts,  the  presence  of  graft-infiltrating 
B  cells  is  a  sign  of  impending  rejection  (25).  Intriguingly,  in 
this  setting,  infiltrating  B  cells  are  uncoupled  from  donor- 
specific  serologic  responses  (29)  and  instead  are  thought  to 
facilitate  T-cell  responses  by  serving  as  APC  and  by  releasing 
cytokines  and  chemokines  that  recruit  T  cells  to  affected 
tissues.  Indeed,  treatment  of  renal  allograft  recipients  with 
anti-CD20  antibody  (rituximab)  not  only  depletes  B  cells 
but  also  results  in  reduced  expression  of  T-cell-associated 
genes  (24).  Similar  mechanisms  are  thought  to  operate  in 
autoimmunity.  In  a  murine  xenograft  model  of  rheumatoid 
arthritis,  depletion  of  B  cells  with  anti-CD20  antibody  lead 
to  the  disappearance  of  not  only  infiltrating  B  cells  but  also 
reduced  activation  and  infiltration  of  T  cells  (27).  These 
findings  raise  the  intriguing  possibility  that  CD20+  TIL  may 
engage  in  similar  cooperative  interactions  with  T  cells  in 
cancer. 

In  this  study,  we  investigated  the  possible  functional  roles 
of  CD20+  TIL  in  HGSC.  Using  matched  tumor  and  serum 
samples  from  HGSC  patients,  we  examined  the  relationship 
between  CD20+  TIL  and  tumor-specific  autoantibody 
responses.  Using  multiparameter  immunohistochemistry 


(IHC)  and  flow  cytometry,  we  assessed  potential  contribu¬ 
tions  of  CD20+  TIL  to  cellular  tumor  immunity.  Our  results 
provide  evidence  that  CD20+  and  CD8+  TIL  work  cooper¬ 
atively  to  mediate  antitumor  immunity  in  HGSC,  leading  to 
markedly  prolonged  patient  survival. 

Materials  and  Methods 

Study  subjects  and  sample  processing 

Studies  were  carried  out  with  3  HGSC  patient  cohorts. 
Cohort  A  (n  =  40),  assembled  through  a  prospective  study 
entitled  Immune  Response  to  Ovarian  Cancer  (IROC), 
included  serial  serum  samples,  cryopreserved  cells,  frozen 
tissue,  and  formalin-fixed  paraffin-embedded  (FFPE)  tissue 
blocks.  Clinical  characteristics  are  described  in  Supplemen¬ 
tary  Table  SI .  Cohort  B  was  a  small  retrospective  cohort  (n  = 
30)  with  matched  serum  samples  and  FFPE  tissue  blocks. 
The  clinical  characteristics  were  similar  to  those  of  cohort  A. 
Cohort  C  (n  =  194)  was  from  a  large  previously  published 
retrospective  cohort,  which  included  only  FFPE  tissue 
blocks  (6,  30).  The  clinical  characteristics  were  similar  to 
those  of  cohort  A,  with  the  exception  that  all  patients  were 
optimally  debulked  during  primary  surgery.  For  ELISA 
experiments,  the  control  group  consisted  of  age-matched 
women  with  both  ovaries  intact  and  no  personal  history  of 
cancer.  All  specimens  and  clinical  data  were  collected  with 
informed  patient  consent  through  the  BC  Cancer  Agency's 
Tumour  Tissue  Repository  (cohort  A)  or  OvCaRe  (cohorts  B 
and  C)  under  protocols  approved  by  the  Research  Ethics 
Board  of  the  BC  Cancer  Agency  and  the  University  of  British 
Columbia. 

Tumor  tissue  was  obtained  at  the  time  of  primary  cytor¬ 
eductive  surgery  before  any  other  treatment.  Samples  were 
either  preserved  in  formalin  and  processed  in  paraffin  for 
preparation  of  tissue  microarrays,  frozen  in  OCT  medium 
for  sectioning  and  DNA  extraction,  or  processed  into  single¬ 
cell  suspensions  as  previously  described  (31).  Serum  was 
processed  using  serum  separator  tubes,  aliquoted,  and 
stored  at  —  80° C.  Pretreatment  serum  was  collected  before 
primary  surgery  (except  for  patient  IROC037,  which  was 
collected  1  month  postsurgery),  and  subsequent  samples 
were  collected  approximately  3,  6,  12,  and  24  months 
following  surgery.  Peripheral  blood  mononuclear  cells 
(PBMC)  from  patients  and  healthy  controls  were  collected 
in  heparinized  Vacutainer  tubes,  isolated  by  Ficoll  density 
centrifugation  and  cryopreserved. 

Immunohistochemistry 

For  cohort  A,  a  tissue  microarray  (TMA)  was  constructed 
using  a  manual  tissue  arrayer  from  Beecher  Instruments. 
The  TMA  consisted  of  duplicate  1-mm  cores  taken  from 
representative  central  regions  of  nonnecrotic  tumor  epithe¬ 
lium  with  some  stromal  regions,  as  defined  upon  review  of 
hematoxylin  and  eosin-stained  sections  by  a  pathologist 
(PHW).  For  cohorts  B  and  C,  TMAs  were  constructed  in  a 
similar  manner  but  consisted  of  0.6-mm  duplicate  cores 
and  were  constructed  using  a  Pathology  Devices  tissue 
arrayer.  TMAs  from  cohorts  A  and  B  were  prepared  and 
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stained  as  previously  described  (31)  using  an  anti-CD20 
rabbit  polyclonal  primary  antibody  (catalog  #RB-9013;  Lab 
Vision).  The  cohort  C  TMA  was  prepared  and  stained  as 
previously  described  (30). 

Scoring  of  TMAs  for  cohorts  A  and  B  was  carried  out  in  a 
blinded  fashion  by  a  pathologist  (PHW).  Tumor  cores  were 
initially  scored  by  assessing  at  low  magnification  the  pro¬ 
portion  of  each  core  that  comprised  tumor  or  stroma  and 
then  scoring  at  20  x  magnification  the  number  of  positively 
stained  TIL  within  the  core  area  (through  direct  counting  up 
to  20  or  by  estimation  in  10s  when  in  excess  of  this 
number),  as  well  as  the  proportion  of  TIL  present  in  tumor 
epithelium  versus  stroma.  For  cohort  C,  cores  were  scored 
either  0  (no  cells  present),  1(1-5  cells),  2  (6-19  cells),  or  3 
(20+  cells)  and  then  binarized  to  negative  versus  positive 
(score  0  vs.  scores  1,  2,  and  3).  This  binarization  was  used  to 
generate  the  survival  data  in  Fig.  6E.  Survival  curves  were  not 
determined  for  cohorts  A  and  B  due  to  small  sample  size  and 
insufficient  follow-up  time. 

Two-color  IHC  was  carried  out  on  FFPE  sections  follow¬ 
ing  deparaffinization  in  a  decloaking  chamber  (Biocare 
Medical)  with  Diva  Decloaker.  Samples  were  blocked  with 
Peroxidased-1  and  Background  Sniper  and  then  labeled 
with  a  cocktail  of  antibodies  to  CD 8  (mouse  monoclonal 
C8/144B;  Cell  Marque)  and  CD20  (rabbit  polyclonal  cat¬ 
alog  #  E2560;  Spring  Bio)  for  30  minutes  at  room  temper¬ 
ature.  The  two-color  polymer  kit  Mach-2  DS2  was  used  to 
amplify  the  primary  signal,  and  Betazoid  DAB  and  Warp 
Red  chromogens  were  used  for  detection.  Slides  were  coun- 
terstained  with  hematoxylin.  All  staining  reagents  were 
from  Biocare  Medical.  Images  were  captured  using  an 
Olympus  BX53  microscope  and  Nuance  multispectral 
imaging  system  (CRI). 

Flow  cytometry 

Cryopreserved  cells  were  labeled  with  antibodies  listed  in 
Supplementary  Table  S2.  Data  were  acquired  using  a  BD 
Influx  cell  sorter  (BD  Biosciences)  and  analyzed  using 
Flowjo  software  (Tree  Star).  All  samples  are  gated  on  lym¬ 
phocytes  based  on  forward  and  side  scatter  plots.  Antibody 
specificity  was  tested  using  matched  isotypes  and  fluores- 
cence-minus-one  controls. 

Idiotype  sequencing 

Six  independent  PCRs  were  carried  out  using  primer  pairs 
designed  to  amplify  VH1,  VH2,  VH3,  VH4,  VH5,  and  VH6 
sequences.  After  amplification,  PCRs  were  combined  in 
pools  of  two  and  cloned  into  the  pENTR/D-TOPO  plasmid 
vector  (Invitrogen).  At  the  BC  Cancer  Agency's  Genome 
Sciences  Centre  (Vancouver,  Canada),  plasmids  were  elec¬ 
troporated  into  DH10B  T1  phage-resistant  cells,  plated 
onto  selective  growth  media,  and  robotically  arrayed  in 
384-well  plates.  Plasmid  DNAs  were  then  extracted  and 
BigDye  Terminator  (ABI)  cycle  sequenced  on  ABI  3730 
sequencers  using  conventional  procedures. 

DNA  sequences  were  analyzed  using  the  IMGT/V-QUEST 
tool  (32)  of  IMGT,  the  international  ImMunoGeneTics 
information  system  (33).  Up  to  96  sequences  were  analyzed 


from  each  tumor  sample.  Sequences  were  compared  with 
germline  sequences  to  identify  somatic  hypermutations, 
and  V-D-J  junctions  were  compared  to  identify  clonal 
sequences. 

ELISA  to  detect  serum  antibodies  to  NY-ESO-1  and  p53 

cDNAs  encoding  the  autoantigens  NY-ESO-1  and  p53 
were  amplified  by  reverse  transcriptase  PCR  (RT-PCR)  from 
the  ovarian  cancer  cell  line  OVCAR3,  cloned  into  the 
pDEST17  vector,  and  used  for  recombinant  protein  pro¬ 
duction,  as  previously  described  (NY-ESO-1;  ref.  31)  or 
using  the  primer  pair:  5'-CACCATGGAGGAGCCGCAGTCA 
and  S'-TTAGTCTGAGTCAGGCCCTrCTGT  (p53).  Serum 
was  assayed  in  triplicate  for  IgG  antibodies  to  recombinant 
NY-ESO-1  and  p53  by  ELISA  as  previously  described  (31). 
Triplicate  measurements  of  control  sera  were  used  to  estab¬ 
lish  the  mean  and  SD  of  autoantibody  responses  to  each 
antigen.  For  cancer  patients,  autoantibody  responses  were 
scored  as  positive  if  they  were  significantly  greater  than  the 
mean  of  the  controls.  Specifically,  the  mean  of  the  controls 
was  subtracted  from  individual  patient  responses,  which 
were  then  divided  by  the  SD  of  all  controls.  Responses 
greater  than  3.0  SDs  from  the  control  mean  were  considered 
positive  (34).  Association  between  CD20  status  and  auto¬ 
antibody  responses  was  assessed  using  the  Fisher  exact  test 
with  a  2-tailed  Rvalue  (GraphPad  Prism).  To  assess  quan¬ 
titative  changes  over  time,  autoantibody  titers  were  deter¬ 
mined  using  a  method  similar  to  Gnjatic  and  colleagues 
(19).  Each  serum  sample  was  subjected  to  4-fold  serial 
dilutions  from  1:100  to  1:6,553,600.  Reciprocal  titers  were 
determined  based  on  the  highest  dilution  that  gave  a 
positive  response  against  NY-ESO-1.  The  cutoff  for  positive 
responses  was  based  on  optical  density  values  greater  than 
twice  the  value  of  the  average  optical  density  of  the  4  lowest 
dilutions  of  a  pool  of  5  healthy  donor  serum  samples. 

Results 

CD20+  TIL  are  mature,  antigen-experienced  B  cells 

We  previously  reported  that  CD20+  TIL  are  strongly 
associated  with  prolonged  survival  in  HGSC  (6).  To  inves¬ 
tigate  potential  immunologic  mechanisms  underlying  this 
association,  we  first  assessed  whether  CD20+  TIL  exhibit 
characteristics  of  antigen-experienced  B  cells.  We  selected  9 
HGSC  tumors  containing  CD20+  TIL,  as  determined  by  IHC 
(Fig.  1  A)  and  used  multiparameter  flow  cytometry  to  assess 
surface  immunoglobulin  expression  on  CD20+  B  cells  from 
tumor  tissue  and  matched  blood  samples.  In  peripheral 
blood,  HGSC  patients  showed  a  normal  distribution  of 
naive  (IgD+IgM+IgG_)  and  activated  (IgD_IgM_IgG+) 
CD20+  B  cells,  which  was  similar  to  cancer-free  controls 
(Fig.  IB).  In  contrast,  the  vast  majority  of  CD20+  TIL 
showed  an  activated  phenotype  (IgD-IgM-IgG+),  indicat¬ 
ing  they  had  undergone  Ig  class  switching  in  response  to 
antigen  exposure  (Fig.  IB). 

To  further  investigate  whether  CD20+  TIL  were  antigen 
experienced,  we  sequenced  up  to  96  IgG  heavy  chain  var¬ 
iable  regions  (including  V,  D,  and  J  gene  segments)  from 
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Figure  1.  Surface  immunoglobulin  expression  by  CD20+  TIL.  A, 
immunohistochemical  analysis  of  CD20+  TIL  in  representative  positive 
(left)  and  negative  (right)  HGSC  tumors.  Scale  bar:  50  |im.  B,  flow 
cytometry  data  from  a  representative  patient  showing  that  the  majority  of 
CD20+  TIL  express  surface  IgG,  whereas  CD20+  B  cells  from  peripheral 
blood  (PBMC)  of  patients  and  healthy  controls  predominantly  express 
IgD  and  IgM.  Samples  were  gated  on  CD20+  lymphocytes  (n  =  9). 

tumor-infiltrating  B  cells  isolated  from  3  HGSC  patients.  We 
observed  a  diverse  cross-section  of  sequences  from  all  major 
VH  families  (Supplementary  Fig.  SI).  Moreover,  tumor- 
infiltrating  B  cells  exhibited  a  high  level  of  somatic  hyper¬ 
mutation:  on  average  20  nucleotides  were  mutated  in  each 
approximately  400  nucleotide  variable  region,  resulting  in 
an  average  of  13  amino  acid  (aa)  substitutions  (Fig.  2A). 
This  is  consistent  with  reference  values  for  antigen-experi¬ 
enced  B  cells  (~18  nucleotide  and  10  aa  substitutions) 
versus  naive  B  cells  (~3  nucleotide  and  3  aa  substitutions; 
ref.  18).  Moreover,  in  all  3  tumors,  there  was  evidence  of  B 
cells  actively  undergoing  somatic  hypermutation,  as  numer¬ 
ous  sequences  were  identified  with  identical  V-D-J  junctions 
but  distinct  mutation  patterns.  An  example  of  a  clone  with 
multiple  mutated  sequences  is  depicted  in  Fig.  2B. 

Finally,  the  clonality  of  CD20+  TIL  was  estimated  by 
calculating  the  prevalence  of  each  variable  gene  sequence  in 
the  dataset.  Sequences  were  defined  as  clonal  if  at  least  2 
independent  sequences  contained  identical  V-D-J  junctions 
but  differing  mutational  patterns.  Using  these  criteria,  we 
found  evidence  for  at  least  11  to  14  B  cell  clones  in  each 
tumor,  some  of  which  were  highly  represented.  For  exam¬ 
ple,  in  IROC015,  16%  of  sequences  were  derived  from  a 
single  B  cell  clone,  an  additional  17%  of  sequences  were 
derived  from  2  other  predominant  clones,  a  further  33%  of 
sequences  were  derived  from  8  minor  clones,  and  the 
remaining  34%  of  sequences  were  unrelated  (Fig.  2C). 
These  trends  were  very  similar  across  all  3  patients  (Fig. 
2C),  indicating  that  CD20+  TIL  represent  oligoclonal  popu¬ 
lations.  Note  that  these  are  conservative  estimates  of  clonal 


abundance,  as  we  disregarded  any  identical  sequences  (i.e., 
sequences  that  had  the  same  mutational  profiles),  reason¬ 
ing  that  these  could  result  from  RT-PCR  amplification  of 
mRNA  molecules  derived  from  the  same  B  cell.  When  the 
analysis  was  repeated  with  the  inclusion  of  duplicate 
sequences,  the  predominant  B-cell  clones  in  each  patient 
represented  22%  to  24%  of  all  sequences  (compared  with 
8%-16%  of  sequences  in  the  initial  analysis).  We  saw  no 
evidence  of  shared  clonotypes  between  patients  (Supple¬ 
mentary  Fig.  SI).  Collectively,  these  results  indicated  that 
CD20+  TIL  have  undergone  activation,  Ig  class  switching, 
somatic  hypermutation,  and  clonal  expansion,  all  of  which 
are  hallmarks  of  antigen  exposure. 

CD20+  TIL  have  an  atypical  CD2  7-negative  memory 
phenotype 

Given  the  above  results,  we  further  characterized  the 
phenotype  of  antigen-experienced  CD20+  TIL  by  flow  cyto¬ 
metry  to  determine  whether  they  exhibit  features  of  germi¬ 
nal  center  or  memory  B  cells.  Triple  staining  with  antibodies 
against  CD20,  CD38,  and  IgD  revealed  that  CD20+ 
PBMC  from  HGSC  patients  and  healthy  controls  consisted 
of  2  populations:  50%  to  65%  of  CD20+  PBMC  were 
IgD+CD38“/low,  indicative  of  naive  B  cells,  whereas 
35%  to  50%  were  IgD  “CD  3  8  “/low,  indicative  of  memory 
B  cells  (Fig.  3A,  top  panels).  The  results  for  CD20+  TIL 
were  strikingly  different.  Naive  CD20+  B  cells  were  not 
seen.  In  3  of  9  TIL  samples,  a  small  population  of 
IgD  “CD  3  8+  cells  was  observed,  consistent  with  germinal 
center  B  cells  (Fig.  3A,  rightmost  CD20+  TIL  panel).  Most 
notably,  in  all  9  patients  the  vast  majority  of  CD20+  TIL 
were  IgD“CD38“^low,  indicative  of  a  memory  B-cell  phe¬ 
notype  (Fig.  3A). 

CD27  is  a  cell-surface  protein  that  is  expressed  under 
normal  physiologic  conditions  (35-37).  As  expected,  the 
vast  majority  of  peripheral  blood  IgG+  B  cells  from  HGSC 
patients  and  healthy  controls  expressed  CD27  (Fig.  3B), 
consistent  with  a  memory  phenotype,  whereas  IgD+  (naive) 
B  cells  lacked  CD27  expression  (data  not  shown).  In  sharp 
contrast,  in  8  of  9  tumor  samples,  virtually  all  IgG+CD20+ 
TIL  failed  to  express  CD27  (Fig.  3B),  suggesting  that  they 
represent  atypical  CD27“  memory  B  cells,  as  described  in 
other  physiologic  settings  (38-41). 

In  addition  to  CD20+  B  cells,  we  analyzed  tumor 
and  blood  samples  for  the  presence  of  plasma  cells 
and  plasmablasts.  There  was  a  complete  absence  of 
CD20“CD38+CD138+  plasma  cells  in  all  analyzed  tumors, 
whereas  the  majority  of  analyzed  tumors  (80%)  contained 
small  numbers  of  cells  with  a  plasmablast  phenotype 
(CD20“CD38+CD138“CD3“CD56“;  Fig.  3C).  In  summa¬ 
ry,  B-cell  lineage  TIL  in  HGSC  are  composed  predominantly 
of  CD27“  memory  B  cells,  with  small  numbers  of  germinal 
center  B  cells  and  plasmablasts. 

CD20+  TIL  are  uncoupled  from  tumor- specific  serum 
autoantibody  responses 

The  low  abundance  of  plasmablasts  and  plasma  cells  in 
tumors  suggested  that  CD20+  TIL,  though  antigen- 
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Figure  2.  CD20+ TIL  are  oligoclonal 
and  have  undergone  somatic 
hypermutation.  A,  by  DNA  sequence 
analysis  IgG  heavy  chain  genes  from 
3  HGSC  tumors  displayed  extensive 
somatic  hypermutation.  Solid  lines: 
expected  values  for  antigen- 
experiened  B  cells,  dotted  lines: 
expected  values  for  naive  B  cells  (1 8). 
B,  a  representative  IgG  clonal  group 
showing  clonal  evolution  within  the 
tumor.  The  germline  sequence  is 
displayed  on  the  bottom  row, 
whereas  all  other  rows  represent 
unique  sequences  present  in  IgG 
molecules  derived  from  tumor  tissue. 
Black  diamonds:  nonsynonymous 
mutations,  white  diamonds: 
synonymous  mutations,  letters: 
amino  acid  substitutions.  C,  IgG 
sequences  from  3  HGSC  tumors 
reveal  the  oligoclonal  composition  of 
CD20+  TIL.  Eight  to  16  percent  of 
sequences  were  derived  from  a 
single  "primary"  clone  in  each  tumor, 
whereas  a  large  percentage  of  the 
remaining  sequences  were  derived 
from  other  predominant  or  minor 
clones. 
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experienced,  might  not  contribute  significantly  to  systemic 
humoral  immunity.  To  further  explore  this  possibility,  we 
evaluated  whether  CD20+  TIL  were  correlated  with  tumor- 
specific  serum  autoantibody  responses.  We  and  others  have 
shown  that  at  least  40%  of  HGSC  patients  have  serum 
autoantibody  responses  to  one  or  more  tumor  antigens 
(19,  20).  Of  the  target  antigens  identified  to  date,  NY- 
ESO-1  and  p53  are  among  the  most  frequently  recognized 
in  ovarian  cancer  (19,  20),  so  these  2  antigens  were  selected 
for  further  analysis.  By  ELISA,  8  of  40  (20%)  HGSC  patients 
in  cohort  A  showed  a  serum  autoantibody  response  to  NY- 


ESO-1  or  p53  (Fig.  4Aand  B).  Notably,  these  responses  were 
equally  prevalent  in  patients  with  or  without  CD20+  TIL 
(19%  vs.  20%,  respectively;  P  =  1.00;  Fig.  4Aand  B).  Similar 
results  were  seen  in  an  independent  cohort  of  HGSC 
patients  (n  =  30;  cohort  B),  in  which  autoantibody 
responses  were  equally  prevalent  in  patients  with  or  without 
CD20+  TIL  (30%  vs.  29%,  respectively;  P  =  1.00;  data  not 
shown). 

To  further  address  this  issue,  we  assessed  whether 
removal  of  CD20+  TIL,  which  occurs  during  primary  sur¬ 
gery  and  chemotherapy,  abrogated  tumor-specific  serum 
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Figure  3.  Flow  cytometric  characterization  of  B  cells,  plasma  cells,  and  plasmablasts  in  HGSC.  A,  naive  versus  memory  B  cells:  Flow  cytometric  analysis  of 
PBMC  and  TIL  samples  stained  with  antibodies  to  CD20,  IgD,  and  CD38.  As  expected,  CD20+  PBMC  from  healthy  controls  and  HGSC  patients  consist 
primarily  of  naive  (lgD+CD38~)  and  memory  (IgD  CD38-)  B  cells  (top).  In  contrast,  the  vast  majority  of  CD20+  TIL  show  a  memory  phenotype  (IgD  CD38  ; 
bottom).  Data  are  representative  of  9  HGSC  patients.  In  3  of  9  patients  (as  shown  for  IROC015),  a  small  population  of  germinal  center  CD20+  B  cells 
(IgD  CD38+)  is  also  seen.  B,  CD27  expression:  Flow  cytometric  analysis  of  CD27  expression  by  PBMC  and  TILfrom  a  representative  healthy  donor  and  HGSC 
patient.  Samples  are  gated  on  CD20+lgG+  lymphocytes.  As  expected,  the  majority  of  CD20+lgG+  B  cells  in  PBMC  express  CD27.  In  contrast, 
the  vast  majority  of  CD20+lgG+  TIL  fail  to  express  CD27.  Data  are  representative  of  9  HGSC  patients.  C,  Plasmablasts  and  plasma  cells:  Flow 
cytometric  analysis  of  CD38  and  CD138  expression  by  PBMC  and  TIL  from  a  representative  healthy  donor  and  a  HGSC  patient  positive  for  CD20+  TIL. 
Samples  were  gated  on  CD45+CD20  lymphocytes.  Left,  PBMC  and  TIL  samples  generally  lack  CD20  CD38+CD138+  plasma  cells  (expected  to  be  in  top 
right  quadrant).  Right,  further  gating  of  CD38+  cells  from  the  left  panels  reveals  the  presence  of  low  numbers  of  cells  with  a  plasmablast  phenotype 
(CD20-CD38+CD3-CD56  ,  bottom  left  quadrants)  in  PBMC  and  TIL  samples.  Data  are  representative  of  4  of  5  tumor  samples. 


autoantibody  responses.  We  quantified  autoantibody  titers 
to  NY-ESO-1  in  6  HGSC  patients  using  serum  samples 
collected  before  and  2  to  6  months  after  primary  surgery. 
We  selected  serum  samples  that  were  collected  as  near  as 
possible  to  the  time  of  lowest  tumor  burden  according  to 
CA125  levels  (Supplementary  Table  S3).  Five  of  6  patients 
showed  an  objective  response  to  surgery  and  chemothera¬ 
py,  as  manifested  by  CA125  levels,  CT  scan,  or  physical 
examination.  In  most  patients,  auto  antibody  levels 


decreased  after  treatment,  although  one  patient  (IROC065) 
showed  a  reproducible  increase  (Fig.  4C).  Irrespective  of 
these  changes,  all  patients  maintained  high-titer  NY-ESO- 
1 -specific  autoantibody  responses  (range:  1:3,200  to 
1:204,800;  Supplementary  Table  S3)  after  cytoreductive 
treatment.  This  indicated  that  the  majority  of  anti-NY- 
ESO-1  autoantibody  production  occurs  outside  of  tumor 
tissue.  Thus,  by  multiple  lines  of  evidence,  CD20+  TIE  are 
uncoupled  from  humoral  immune  responses  in  HGSC. 
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Figure  4.  CD20+  TIL  are  uncoupled  from  humoral  immunity.  A  and  B, 
ELISA  data  from  cohort  A  (n  =  40)  showing  equivalent  frequencies  of 
autoantibody  responses  to  (A)  NY-ESO-1  and  (B)  p53  in  pretreatment 
sera  from  patients  with  tumors  scored  as  positive  (black  bars)  or  negative 
(white  bars)  for  CD20+  TIL.  The  threshold  is  set  at  3.0  SDs  from  the  mean 
optical  density  of  30  healthy  control  samples.  C,  ELISA  data  for  6  HGSC 
patients  showing  serum  autoantibody  responses  to  NY-ESO-1  before 
primary  treatment  (Pre;  white  bars)  and  2  to  6  months  after  primary 
surgery  when  tumor  burden  was  lowest  according  to  CA-1 25  levels  (Post; 
black  bars).  Data  are  plotted  as  reciprocal  titers.  Representative  of 
3  independent  experiments,  each  run  in  triplicate. 

CD20+  TIL  express  cell  surface  markers  characteristic 
of  antigen-presenting  cells 

The  foregoing  results  suggested  that  CD20+  TIL  might 
contribute  to  tumor  immunity  through  cellular  rather  than 
humoral  mechanisms.  Indeed,  in  addition  to  producing 


antibodies,  B  cells  can  facilitate  T-cell  responses  by  serving 
as  APC  and  by  recruiting  T  cells  to  affected  tissues  (42).  We 
therefore  assessed  whether  CD20+  TIL  might  serve  as  APC  to 
T  cells,  as  described  in  other  settings  (23,  26,  42,  43). 
Unfortunately,  the  low  absolute  numbers  of  CD20+  TIL  in 
human  specimens  precluded  in  vitro  functional  analyses. 
However,  sufficient  material  was  available  to  carry  out 
multiparameter  flow  cytometry  to  determine  whether 
CD20+  TIL  express  cell-surface  molecules  involved  in  anti¬ 
gen  presentation.  CD20+  TIL  expressed  MHC  class  I  (pan- 
HLA-A,  B,  and  C)  and  class  II  (HLA-DR)  molecules,  as  well 
as  CD40  (Fig.  5A).  Furthermore,  CD20+  TIL  expressed  the 
costimulatory  molecules  B7-1  (CD80)  and  B7-2  (CD86), 
albeit  at  moderate  levels  (Fig.  5 A) .  Accordingly,  the  majority 
of  tumor-infiltrating  CD4+  and  CD8+  T  cells  displayed  an 
activated  effector  phenotype  (HLA-DR+CD45RO+CD45- 
RATCD62L-;  Fig.  5B).  Collectively,  our  results  are  consis¬ 
tent  with  the  possibility  that  CD20+  TIL  serve  as  APC  to  T 
cells  in  the  tumor  environment. 

CD8+  and  CD20+  TIL  colocalize  in  HGSC 

If  B  cells  serve  as  APC  in  the  tumor  environment,  they 
would  be  expected  to  be  located  in  close  proximity  to  T  cells. 
To  assess  this,  we  carried  out  2-color  IHC  with  antibodies  to 
CD8  and  CD20.  As  previously  reported,  CD8+  and  CD20+ 
TIL  were  found  in  both  tumor  stroma  and  epithelium  (Fig. 
6A  and  B).  Intriguingly,  CD8+  and  CD20+  lymphocytes 
were  frequently  seen  in  close  proximity,  often  in  loosely 
structured  aggregates  (Fig.  6A).  Such  aggregates  were  even 
found  in  tumors  with  low  densities  of  TIL,  suggesting  the 
colocalization  of  CD20+  and  CD8+  TIL  was  an  active 
process.  There  were  numerous  examples  in  which  CD8+ 
and  CD20+  TIL  were  directly  juxtaposed,  consistent  with 
cell-cell  interactions  (Fig.  6C  and  D).  In  summary,  these 
imaging  studies  indicated  that  CD8+  and  CD20+  TIL  are 
frequently  colocalized  in  HGSC,  further  supporting  a  pos¬ 
sible  role  for  CD20+  TIL  as  APC. 

The  combination  of  CD8+  and  CD20+  TIL  is  associated 
with  prolonged  survival  in  HGSC 

Finally,  we  used  data  from  a  large  previously  published 
cohort  (6)  to  assess  the  relationship  between  CD20+  and 
CD8+  TIL  and  patient  survival.  CD20+  and  CD8+  TIL  were 
found  in  41%  and  84%  of  cases,  respectively.  Of  the  cases 
that  were  positive  for  CD20+  TIL,  almost  all  were  also 
positive  for  CD8+  TIL.  Conversely,  of  the  cases  that  were 
positive  for  CD8+  TIL,  approximately  half  also  contained 
CD20+  TIL.  By  Kaplan-Meier  analysis,  cases  that  were 
positive  for  both  CD8+  and  CD20+  TIL  showed  markedly 
greater  disease-specific  survival  compared  with  those  pos¬ 
itive  for  CD8+  TIL  alone  (Fig.  6E).  Thus,  it  seems  that 
CD20+  TIL  can  potently  enhance  the  antitumor  effect  of 
CD8+  TIL  in  ovarian  cancer. 

Discussion 

We  have  investigated  potential  mechanisms  by  which 
CD20+  TIL  contribute  to  tumor  immunity  in  HGSC, 
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Figure  5.  CD20+  TIL  have  a 
surface  phenotype  characteristic 
of  APC,  whereas  CD4+  and  CD8+ 
TIL  have  a  surface  phenotype 
characteristic  of  activated  T  cells. 
A,  representative  flow  cytometry 
data  showing  expression  of  MHC 
class  I  (pan-HLA-A,  B,  C),  MHC 
class  II  (pan-HLA-DR),  CD40, 
CD80,  and  CD86  in  CD20+  B 
cells  from  healthy  control  PBMC 
(top),  HGSC  patient  PBMC 
(middle),  and  TIL  (bottom). 

Black  lines:  specific  staining 
with  the  indicated  antibody. 

Gray  lines:  negative  controls 
lacking  the  indicated  antibody. 
Samples  were  gated  on 
CD45+CD20+  lymphocytes 
(n  =  9).  B,  representative  flow 
cytometry  data  showing  an 
activated  (HLA-DR+ 

CD45RA  CD45RO+CD62L  ) 
surface  phenotype  on  the 
majority  of  CD4+  (left)  and  CD8+ 
(right)  TIL.  Samples  are  gated  on 
CD3+  lymphocytes  (n  =  9). 


thereby  promoting  prolonged  patient  survival.  CD20+  TIL 
showed  clear  evidence  of  antigen  exposure,  including 
expression  of  activation  markers,  class  switching  to  IgG, 
somatic  hypermutation,  and  oligoclonality.  The  majority 
failed  to  express  CD27,  suggesting  that  they  belong  to  a 
recently  described  subset  of  CD27~  memory  B  cells  (39- 
41,  44).  Despite  expressing  IgG,  CD20+  TIL  showed  no 
correlation  with  tumor-specific  serum  autoantibodies  and 
hence  seemed  to  be  uncoupled  from  humoral  immunity. 
Instead,  CD20+  TIL  correlated  strongly  with  CD8+  TIL.  The 
2  lymphocyte  subsets  were  found  in  close  proximity  in 
tumor  tissue,  often  in  loose  aggregates.  CD20+  TIL 
expressed  surface  markers  characteristic  of  APC,  including 
MHC  class  I,  MHC  class  II,  CD80,  and  CD86.  Furthermore, 
tumors  positive  for  both  CD8+  and  CD20+  TIL  were  asso¬ 


ciated  with  prolonged  survival  compared  with  those  posi¬ 
tive  for  CD8+  TIL  alone.  Collectively,  our  results  support  the 
novel  concept  that  CD20+  TIL  enhance  tumor  immunity 
through  functional  interactions  with  CD8+  TIL. 

Sequencing  of  IgG  molecules  from  CD20+  TIL  yielded 
clear  evidence  of  oligoclonality  and  somatic  hypermuta¬ 
tion,  which  are  hallmarks  of  antigen  exposure.  The  degree  of 
B-cell  clonality  we  measured  in  HGSC  is  similar  to  prior 
studies  in  breast  cancer  and  germ  cell  tumors.  Specifically, 
we  found  evidence  for  11  to  14  B-cell  clones  in  HGSC, 
compared  with  6  to  1 3  clones  in  germ  cell  tumors  (18)  and  3 
to  6  clones  in  breast  cancer  (12,  17).  Similarly,  58%  to  66% 
of  sequences  were  clonally  derived  in  HGSC,  compared 
with  18%  to  79%  in  germ  cell  tumors  and  30%  to  69%  in 
breast  cancer.  Notably,  these  are  minimal  estimates,  as  it  is 
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Figure  6.  CD20+  and  CD8+  TIL 
colocalize  and  are  associated  with 
prolonged  survival  in  HGSC.  A,  IHC 
of  HGSC  tumor  (IROC006)  showing 
CD20+  TIL  (red)  and  CD8+  TIL 
(brown)  primarily  in  stromal  regions. 
B,  IHC  of  a  densely  infiltrated  HGSC 
tumor  (IROC042)  showing  CD20+  TIL 
(red)  and  CD8+  TIL  (brown)  in 
epithelial  and  stromal  regions.  C, 
high  power  IHC  image  of  a  HGSC 
tumor  (I ROC01 0)  showing  CD20+  TIL 
(red)  and  CD8+  TIL  (brown)  in  close 
proximity.  D,  the  same  image  as 
shown  in  (C),  but  with  false  coloring 
to  further  contrast  the  CD20  (red)  and 
CD8  (green)  signals.  All  samples 
were  counterstained  with 
hemotoxylin  (blue/purple).  Original 
images  were  taken  with  a  20  x  (A  and 
B)  or  40  x  objective  (C  and  D).  Scale 
bars:  1 00  (im  (A  and  B)  or  50  |im 
(C  and  D).  Representative  of  5  tumor 
samples.  E,  Kaplan-Meier  curves 
showing  that  the  presence  of  both 
CD8+  and  CD20+  TIL  is  associated 
with  prolonged  disease-specific 
survival  compared  with  CD8+  TIL 
alone  or  lack  of  CD8+  TIL.  Primary 
data  is  from  cohort  C  (6).  Log-rank 
tests  were  used  to  derive  P  values  for 
comparisons  between  each  group 
and  between  all  3  groups. 
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possible  that  additional  less  abundant  clones  were  present 
in  all  of  these  studies.  In  addition  to  clonality,  the  degree  of 
somatic  hypermutation  we  observed  in  HGSC  is  similar  to 
that  reported  for  breast  cancer  and  germ  cell  tumors. 

What  antigens  might  be  recognized  by  CD20+  TIL? 
Although  this  issue  has  yet  to  be  addressed  in  HGSC,  prior 
studies  in  breast  cancer  identified  2  antigens  recognized  by 
recombinant  antibodies  derived  from  tumor-infiltrating  B 
cells:  ganglioside  D3  and  a  cleaved  form  of  (3-actin  that  was 
shown  to  become  exposed  on  the  surface  of  apoptotic  cells 
(13,  45).  These  results  suggest  that  CD20+  TIL  responses 
reflect  a  breakdown  of  tolerance  to  self-proteins,  possibly 
due  to  apoptosis  of  tumor  cells.  Identification  of  additional 
antigens  recognized  by  CD20+  TIL  may  yield  further 


insights  into  the  mechanisms  by  which  B-cell  responses  are 
triggered  and  contribute  to  tumor  immunity. 

Even  though  CD20+  TIL  express  surface  immunoglobu¬ 
lin,  we  found  that  they  are  not  associated  with  serum 
autoantibodies  to  the  common  tumor  antigens  p53  and 
NY-ESO-1.  This  lack  of  correlation  was  seen  across  2  inde¬ 
pendent  cohorts  of  patients,  as  well  as  within  individual 
patients  undergoing  cytoreductive  surgery  and  chemother¬ 
apy.  Although  only  2  antigens  were  assessed,  these  are  the 
most  frequently  recognized  tumor  antigens  in  HGSC  iden¬ 
tified  to  date  ( 1 9,  20) .  This  finding  is  consistent  with  the  fact 
that  CD20+  TIL  are  associated  with  prolonged  survival  in 
HGSC,  whereas  serum  autoantibodies  show  little  or  no 
positive  prognostic  signal  in  most  studies  (31,  46-48). 
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Collectively,  our  data  supports  the  conclusion  that  serum 
autoantibodies  are  produced  by  cells  outside  the  tumor. 
Prime  candidates  are  the  plasma  cells  in  bone  marrow, 
which  are  the  primary  source  of  serum  antibodies  in  healthy 
individuals.  It  remains  possible  that  CD20+  TIL,  and  the 
small  populations  of  germinal  center-like  B  cells  and 
CD38+  plasmablasts  we  identified  in  tumors,  are  pre¬ 
cursors  to  plasma  cells  that  in  turn  produce  serum  anti¬ 
bodies  in  some  patients. 

Although  CD20+  TIL  exhibit  the  hallmarks  of  mature, 
antigen-experienced  B  cells,  the  majority  failed  to  express 
the  conventional  memory  B-cell  marker  CD27.  CD27- 
memory  B  cells  are  present  at  low  numbers  in  healthy 
donors  (39-41,  44).  In  patients  with  systemic  lupus  erythe¬ 
matosus,  the  number  of  CD27-  memory  B  cells  is  elevated 
in  proportion  to  disease  severity  (44).  To  our  knowledge, 
this  study  is  the  first  to  identify  CD27-  memory  B  cells  in 
cancer.  Although  the  functional  differences  between  CD27+ 
and  CD27-  memory  B  cells  have  yet  to  be  elucidated, 
several  possibilities  can  be  imagined.  Activated  helper  T 
cells  express  CD 70,  which  can  ligate  CD27  on  B  cells  to 
induce  Ig  secretion  (35).  Perhaps  the  lack  of  CD 2 7  expres¬ 
sion  on  CD20+  TIL  allows  them  to  forego  antibody 
production  and  instead  contribute  to  cellular  immunity. 
Alternatively,  the  lack  of  CD 2 7  expression  could  indicate 
recent  activation  of  B  cells,  as  exposure  of  B  cells  to  CD70+T 
cells  leads  to  rapid  CD27  downregulation  (35).  Perhaps  the 
presence  of  CD27-  memory  B  cells  reflects  a  response  to 
CD70  signals  by  CD4+  and/or  CD8+  TIL  in  HGSC. 

Using  two-color  IHC,  we  showed  that  CD20+  TIL  are 
often  found  in  close  proximity  to  CD8+  TIL  (Fig.  6A-D). 
Moreover,  the  presence  of  both  CD20+  and  CD8+  TIL  was 
associated  with  markedly  prolonged  patient  survival  com¬ 
pared  with  CD8+  TIL  alone  (Fig.  6E).  These  findings  suggest 
that  CD20+  TIL  may  promote  cytolytic  antitumor 
responses.  In  theory,  this  could  occur  through  multiple 
mechanisms.  First,  CD20+  TIL  could  help  recruit  and  retain 
T  cells  at  the  tumor  site,  facilitating  the  formation  of  the 
lymphoid  aggregates  we  observed  (Fig.  6A-D).  Indeed,  in 
autoimmunity,  allograft  rejection,  and  chronic  infection,  B 
cells  are  central  to  the  formation  of  so-called  tertiary  lym¬ 
phoid  structures  in  affected  tissues  (28,  49).  For  example,  in 
a  xenograft  model  of  rheumatoid  arthritis,  depletion  of  B 
cells  with  anti-CD20  antibody  led  to  decreased  T-cell  acti¬ 
vation  and  infiltration  (27).  Conversely,  depletion  ofT  cells 
led  to  disruption  of  tertiary  lymphoid  structures  and  loss  of 
immunoglobulin  production  by  infiltrating  B  cells  (50). 
Thus,  during  prolonged  immune  responses,  B  cells  and  T 
cells  can  engage  in  cooperative  interactions  to  maintain  a 
strong  presence  at  the  affected  site. 

Second,  CD20+  TIL  could  serve  as  APC  in  the  tumor 
environment.  It  is  well  documented  that  B  cells  can  present 
antigens  to  T  cells,  including  cross-presentation  to  CD8+  T 
cells  (5 1  -53) .  B  cells  have  been  shown  to  serve  as  APC  under 
a  wide  variety  of  physiologic  conditions,  including  auto¬ 
immunity  and  allograft  rejection.  For  example,  B  cells  serve 
as  APC  to  induce  T-cell  responses  in  nonobese  diabetic  mice 
(26).  Similarly,  B  cells  from  multiple  sclerosis  patients  have 


been  shown  to  present  myelin  basic  protein  to  CD4+  T  cells 
in  vitro  (43).  Thus,  B  cells  have  the  capacity  to  function  as 
APC  in  multiple  autoimmune  disorders  and  may  function 
similarly  in  antitumor  responses.  In  support  of  this  notion, 
we  found  that  CD20+  TIL  express  MHC  class  I  and  II,  as  well 
as  the  costimulatory  molecules  CD 80  and  CD 86.  Although 
speculative,  the  presence  of  professional  APC  in  the  tumor 
environment  may  allow  for  sustained  T-cell  responses,  for 
example,  by  providing  supportive  niches  for  memory  T-cell 
proliferation.  Notably,  CD20+  TIL  are  far  more  numerous 
than  immature  or  mature  dendritic  cells  (as  defined  by 
CD  la  or  CD208  expression,  respectively)  in  immune  infil¬ 
trates  in  HGSC  (6).  Moreover,  compared  with  dendritic 
cells,  B  cells  have  the  advantage  that  B-cell  receptor-medi¬ 
ated  endocytosis  facilitates  concentration  of  small  quanti¬ 
ties  of  specific  antigens,  which  may  allow  the  amplification 
of  responses  to  tumor  antigens  expressed  at  low  levels.  To 
definitively  show  that  CD20+  TIL  serve  as  APC  in  the  tumor 
environment,  it  will  be  necessary  to  identify  their  cognate 
antigens  and  assess  recognition  by  CD4+  and  CD8+  TIL. 

Our  findings  have  implications  for  the  immunotherapy 
of  cancer.  Most  cancer  vaccines  are  modeled  on  acute  viral 
infections,  which  are  characterized  by  large  yet  transient  T- 
cell  responses.  By  contrast,  in  autoimmunity  or  allograft 
rejection,  T-cell  responses  are  sustained  over  many  months 
or  years  and  involve  cooperative  interactions  with  coloca¬ 
lized  B  cells.  Thus,  it  may  be  possible  to  generate  more 
potent,  sustained  T-cell  responses  in  the  tumor  environ¬ 
ment  by  promoting  the  infiltration  of  tumor-reactive  B  cells. 
In  support  of  this  idea,  adoptive  transfer  of  activated  B  cells 
has  been  shown  to  induce  tumor-specific  T-cell  immunity  in 
murine  models  (54).  Thus,  improved  understanding  of  the 
functional  phenotype  of  CD20+  TIL,  their  target  antigens, 
and  their  mechanism  of  recruitment  to  target  tissues  may 
facilitate  the  design  of  more  effective  immunotherapies  for 
the  treatment  of  cancer. 
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Supplementary  Figure  1:  B  cells  in  HGSC  tumors  exhibit  oligoclonal  and  somatically 
hypermutated  immunoglobulin  sequences.  Immunoglobulin  genes  were  sequenced  from 
three  HGSC  tumors  and  compared  with  germline  sequences.  Identical  sequences  or 
sequences  containing  only  a  single  difference  were  excluded.  Remaining  sequences  were 
defined  as  clonal  if  V-D-J  junctions  were  identical.  Tables  represent  gene  usage  and  clonal 
groups  from  each  tumor. 
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Supplementary  Table  1.  Clinical  characteristics  of  patient  cohort  A. 


No.  of  patients 


Age  at  surgery,  median  (range) 


Debulking  Status 


Optimal  (no  visible  residual) 


Minimal  Residual  (<1cm  residual  disease) 


Suboptimal  (>1cm  residual  disease) 


Chemotherapy  regimen _ 


Carboplatin/paclitaxel 


Carboplatin  only 


No  chemotherapy 


Stage 


2 


3 


Unknown 


Grade 


3 


Status 


Alive 


Dead 


Follow-up  time  in  yrs,  median  (range) 


40 


66  (44-79) 


2  (5%) 


10(25%) 


28  (70%) 


24  (60%) 


15(37.5%) 


1  (2.5%) 


2  (5%) 


31  (77.5%) 


6  (15%) 


1  (2.5%) 


13  (32.5%) 


27  (67.5%) 


1.9  (0.1-3. 8) 


Supplementary  Table  2.  Antibodies  used  for  flow  cytometry. 


Antigen 

Fluorochrome 

Clone 

Supplier 

CD3 

FITC 

HIT3a 

BD  Biosciences,  Mississauga,  ON,  Canada 

CD3 

PE-Cy7 

SK7 

BD  Biosciences 

CD4 

APC-H7 

RPA-T4 

BD  Biosciences 

CD8 

PerCP-Cy5.5 

RPA-T8 

eBioscience,  San  Diego,  CA 

CD20 

APC-H7 

2H7 

BD  Biosciences 

CD27 

PE 

M-T271 

BD  Biosciences 

CD38 

PerCP-Cy5.5 

HIT2 

BD  Biosciences 

CD40 

FITC 

5C3 

BD  Biosciences 

CD45 

PE-Cy7 

HI30 

eBioscience 

CD45RA 

APC 

H 1 1 00 

BD  Biosciences 

CD45RO 

FITC 

UCHL-1 

eBioscience 

CD56 

APC 

B159 

BD  Biosciences 

CD62L 

PE 

Dreg  56 

BD  Biosciences 

CD80 

PE-Cy5 

L307.4 

BD  Biosciences 

CD86 

PE-Cy5 

2331  (FUN-1) 

BD  Biosciences 

CD138 

PE 

DL-101 

eBioscience 

HLA-ABC 

PE 

W6/32 

eBioscience 

HLA-DR 

FITC 

LN3 

eBioscience 

HLA-DR 

PE 

L243  (G46-6) 

BD  Biosciences 

igD 

FITC 

IA6-2 

BD  Biosciences 

igG 

PE-Cy7 

G18-145 

BD  Biosciences 

igM 

APC 

G20-127 

BD  Biosciences 

Supplementary  Table  3.  CA125  levels  and  NY-ESO-1  autoantibody  titres  before 
(Pre)  and  after  (Post)  primary  surgery  and  chemotherapy  for  the  6  HGSC  patients 
shown  in  Figure  4C.  Titres  represent  mean  values  from  three  independent 


experiments. 


Patient 

CD20-TIL 

CA125 

NY-ESO-1  reciprocal  titre 

status 

Pre 

Post 

Pre 

Post 

IROC008 

CD20- 

2,388 

1,096 

12,800 

3,200 

IROC013 

CD20- 

273 

5 

273,067 

6,400 

IROC033 

CD20- 

6,572 

17 

2,867,200 

19,200 

IROC060 

CD20+ 

4,037 

Not  available 

409,600 

204,800 

IROC065 

CD20+ 

592 

57 

1,167 

4,800 

IROC072 

CD20+ 

184 

6 

409,600 

19,200 
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Abstract 

CD 8+  tumor-infiltrating  lymphocytes  (TIL)  are  associated  with  survival  in  a  variety  of  cancers.  A  second 
subpopulation  of  TIL,  defined  by  forkhead  box  protein  P3  (FoxP3)  expression,  has  been  reported  to  inhibit 
tumor  immunity,  resulting  in  decreased  patient  survival.  On  the  basis  of  this  premise,  several  groups  are 
attempting  to  deplete  FoxP3+ T  cells  to  enhance  tumor  immunity.  However,  recent  studies  have  challenged 
this  paradigm  by  showing  that  FoxP3+  T  cells  exhibit  heterogeneous  phenotypes  and,  in  some  cohorts,  are 
associated  with  favorable  prognosis.  These  discrepant  results  could  arise  from  differences  in  study 
methodologies  or  the  biologic  properties  of  specific  cancer  types.  Here,  we  conduct  the  first  systematic 
review  of  the  prognostic  significance  of  FoxP3+  T  cells  across  nonlymphoid  cancers  (58  studies  from  16 
cancers).  We  assessed  antibody  specificity,  cell-scoring  strategy,  multivariate  modeling,  use  of  single 
compared  with  multiple  markers,  and  tumor  site.  Two  factors  proved  important.  First,  when  FoxP3  was 
combined  with  one  additional  marker,  double-positive  T  cells  were  generally  associated  with  poor 
prognosis.  Second,  tumor  site  had  a  major  influence.  FoxP3+  T  cells  were  associated  with  poor  prognosis 
in  hepatocellular  cancer  and  generally  good  prognosis  in  colorectal  cancer,  whereas  other  cancer  types 
were  inconsistent  or  understudied.  We  conclude  that  FoxP3+  T  cells  have  heterogeneous  properties  that 
can  be  discerned  by  the  use  of  additional  markers.  Furthermore,  the  net  biologic  effects  of  FoxP3+  T  cells 
seem  to  depend  on  the  tumor  site,  perhaps  reflecting  microenvironmental  differences.  Thus,  depletion 
of  FoxP3+ T  cells  might  enhance  tumor  immunity  in  some  patient  groups  but  be  detrimental  in  others.  Clin 
Cancer  Res;  18(11);  1-8.  ©2012  AACR. 


Introduction 

Many  studies  across  a  wide  variety  of  human  cancers  have 
shown  a  clear  association  between  the  presence  of  tumor- 
infiltrating  lymphocytes  (TIL)  and  patient  survival  ( 1  -4) .  To 
further  understand  this  phenomenon,  additional  immune 
markers  have  been  used  to  subdivide  CD3+  T  cells  into 
functional  subsets,  with  special  emphasis  on  cytotoxic  (e.g., 
CD8+,  nucleolysin  TIA-1  isoform  p40+)  and  regulatory 
[e.g.,  CD4+,  interleukin  2  receptor  subunit  alpha  (CD25)+, 
FoxP3+]  phenotypes  (3,  5).  Whereas  TIL-expressing  cyto¬ 
toxic  markers  are  generally  associated  with  favorable  prog¬ 
nosis,  TIL-expressing  regulatory  markers  (referred  to  as 
Tregs)  were  initially  reported  to  correlate  with  poor  prog¬ 
nosis  (5).  This  finding  fit  with  the  general  notion  that  Tregs 
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suppress  adaptive  immune  responses  and  led  many  groups 
to  pursue  strategies  to  deplete  Tregs  from  patients  with 
cancer  as  a  means  to  enhance  tumor  immunity  (6-8). 

In  the  past  decade,  much  effort  has  been  devoted  to 
finding  molecular  markers  that  uniquely  define  Tregs.  Ini¬ 
tially,  these  cells  were  characterized  as  CD4+  and  CD25hlgh 
(9).  Further  investigation  revealed  that  Tregs  express  and 
functionally  depend  on  the  transcription  factor  forkhead 
box  protein  P3  (FoxP3;  ref.  10).  Indeed,  humans  and  mice 
that  lack  an  intact  FOXP3  gene  suffer  a  severe  autoimmune 
syndrome  known  as  immune  dysregulation/polyendocri- 
nopathy/enteropathy/X-linked  syndrome  in  humans  or  the 
Scurfy  phenotype  in  mice  (10,  11).  Given  its  essential  role  in 
Treg  development  and  function,  FoxP3  became  a  popular 
single  marker  for  Treg  studies  in  cancer.  Intriguingly,  studies 
of  the  prognostic  value  of  FoxP3+ T  cells  have  lead  to  highly 
discrepant  findings.  In  some  studies,  tumor-infiltrating 
FoxP3+  T  cells  have  been  associated  with  poor  prognosis, 
consistent  with  the  initial  hypothesis  that  FoxP3+  Tregs 
inhibit  antitumor  immunity.  In  contrast,  other  studies  have 
found  that  FoxP3+  T  cells  are  associated  with  a  favorable 
prognosis. 

How  can  these  widely  discrepant  prognostic  claims  be 
explained?  On  the  one  hand,  they  could  reflect  technical 
differences  among  studies,  including  the  specific  FoxP3 
antibody  used,  scoring  strategy,  and  statistical  methods. 
Alternatively,  the  differing  claims  could  reflect  biologic 
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Translational  Relevance 

Although  forkhead  box  protein  P3+  (FoxP3+)  T  cells 
are  conventionally  thought  to  suppress  tumor  immuni¬ 
ty,  this  idea  has  been  challenged  by  recent  studies  show¬ 
ing  that,  in  some  patient  cohorts,  tumor-infiltrating 
FoxP3+  T  cells  are  associated  with  favorable  prognosis. 
To  investigate  this  apparent  discrepancy,  we  did  a  com¬ 
prehensive  review  of  the  literature  on  the  prognostic 
significance  of  tumor-infiltrating  FoxP3+  T  cells  in 
human  cancer.  We  conclude  that  FoxP3  is  inadequate 
as  a  single  functional  or  prognostic  marker.  Moreover, 
the  prognostic  significance  of  FoxP3+  T  cells  can  vary 
according  to  tumor  site.  Thus,  the  original  view  that 
FoxP3+  T  cells  invariably  suppress  tumor  immunity  is 
oversimplified.  We  require  better  understanding  of  the 
functional  subtypes  of  FoxP3+  T  cells  and  their  biologic 
properties  in  different  tumor  microenvironments  if  we 
wish  to  rationally  modulate  their  behavior  to  enhance 
tumor  immunity. 


factors.  For  example,  it  is  conceivable  that  FoxP3+  T  cells 
exhibit  conventional  regulatory  (i.e.,  inhibitory)  properties 
in  some  contexts  but  not  others.  Alternatively,  FoxP3+  T 
cells  may  be  consistently  regulatory  in  nature  but  appear  as 
favorable  prognostic  markers  in  some  cancers  because  of 
their  association  with  tumor-infiltrating  CD 8+  T  cells  or 
other  effectors  (12,  13).  Others  have  suggested  that,  in 
colorectal  and  gastric  cancers,  FoxP3+  T  cells  may  inhibit 
tumor-promoting  inflammatory  responses  to  microbes, 
which  could  explain  their  association  with  favorable  out¬ 
comes  in  these  and  similar  contexts  (14).  Finally,  emerging 
evidence  indicates  that  FoxP3  expression  encompasses  a 
heterogeneous  population  of  cells  that  contain  both  regu¬ 
latory  T  cells,  which  produce  cytokines  such  as  TGF-(3l  and 
interleukin  10,  and  nonregulatory  T  cells,  which  may 
express  interferon  gamma  and  interleukin  17  (15-19; 
reviewed  in  ref.  20).  Given  these  various  possibilities,  it 
seems  reasonable  to  question  whether  depletion  of  Tregs 
based  on  FoxP3  expression  is  likely  to  be  beneficial  or 
detrimental  to  patients  with  cancer. 

To  investigate  this  controversy,  we  did  a  comprehensive 
and  critical  review  of  the  literature  on  tumor-infiltrating 
FoxP3+  T  cells  and  prognosis  in  human  cancer.  Articles  for 
review  were  identified  during  a  PubMed  search  using  the 
terms  "FoxP3"  and  "cancer"  and  were  vetted  by  title  and 
abstract  by  one  of  the  authors  (R.J.  deLeeuw).  Several 
selection  criteria  were  applied.  First,  we  excluded  studies 
of  lymphoid  cancers,  because  the  immunologic  nature  of 
these  malignancies  makes  it  difficult  to  assess  whether 
FoxP3+ T  cells  are  acting  directly  on  tumor  cells  or  indirectly 
on  antitumor  effector  lymphocytes.  Second,  we  excluded 
studies  that  only  correlated  FoxP3+  T  cells  with  late-stage 
disease  as  opposed  to  patient  survival.  Third,  we  included 
only  those  studies  that  measured  FoxP3  expression  by 
immunohistochemistry  (IHC)  or  immunofluorescence  to 


ensure  that  the  intratumoral  location  of  FoxP3+  cells  was 
known.  Finally,  we  reviewed  a  given  data  set  only  once, 
excluding  secondary  or  tertiary  studies  that  referred  to  a 
previously  published  data  set. 

In  the  end,  we  reviewed  58  studies  encompassing  16 
different  cancer  types  (Table  1),  including  bladder  (19), 
breast  (21-28),  cervical  (29,  30),  colorectal  (12,  31-39), 
endometrial  (40-42),  gastric  (14,  43-46),  head  and  neck 
(47),  hepatocellular  (48-52),  lung  (53,  54),  melanoma 
(55-58),  mesothelial  (59),  oral  (4,  60-63),  ovarian  (2, 
3,  13,  64-67),  pancreatic  (68),  renal  (69,  70),  and  vulvar 
cancers  (Supplementary  Table  SI;  ref.  71).  The  reported 
prognostic  value  of  FoxP3+  T  cells  in  these  58  studies 
ranged  from  poor  (n  =  23),  to  neutral  (n  =  23),  to  good 
(n  =  12).  To  better  understand  why  the  prognostic  value  of 
FoxP3+  T  cells  varies  so  widely,  we  assessed  each  study  for 
technical  factors  (including  the  specific  FoxP3  antibody 
used,  scoring  strategy,  and  the  use  of  multivariate  model¬ 
ing)  and  biologic  factors  (including  the  use  of  additional 
markers  to  define  Tregs  and  the  tumor  site  studied). 

Antibody  Specificity 

Different  FoxP3  antibodies  can  yield  different  staining 
patterns,  indicating  that  some  antibodies  may  have  subop- 
timal  sensitivity  or  specificity  (72,  73).  Although  18  of  the 
58  reviewed  studies  failed  to  state  which  specific  FoxP3 
antibody  was  used,  within  the  remaining  40  studies,  1 1 
different  FoxP3  antibodies  were  used  (Table  1).  The  most 
commonly  used  antibody  was  a  monoclonal  designated 
236A/E7.  In  the  23  studies  that  used  236A/E7,  the  prog¬ 
nostic  significance  of  FoxP3+  T  cells  ranged  from  poor  (n  = 
10),  to  neutral  (n  =  8),  to  good  (n  =  5).  Given  that  a  single 
FoxP3  antibody  can  yield  prognostic  results  this  disparate,  it 
seems  that  FoxP3  prognostic  variability  is  not  solely  attrib¬ 
utable  to  antibody  selection. 

Cell-Scoring  Strategy 

We  investigated  4  aspects  of  the  scoring  strategies  used  to 
categorize  tumors  as  positive  or  negative  for  FoxP3+  T  cells: 
cutoff  points,  intratumoral  location,  use  of  tissue  micro¬ 
arrays  compared  with  whole  sections,  and  computerized 
compared  with  manual  counting  (Table  1).  Although  there 
is  no  standard  cutoff  point  for  TIL  studies,  3  2  out  of  5  8  of  the 
reviewed  studies  used  the  median  number  of  FoxP3+ T  cells 
as  the  cutoff  point.  Within  these  32,  the  distribution  among 
poor,  neutral,  and  good  prognostic  claims  was  16, 1 1,  and  5 
studies,  respectively.  The  remaining  studies  used  a  variety  of 
scoring  strategies,  including  the  presence  compared  with 
absence  of  FoxP3+  T  cells,  the  mean  number  of  FoxP3+  T 
cells,  or  other  criteria.  A  fairly  even  distribution  among 
poor,  neutral,  and  good  prognostic  claims  was  observed 
regardless  of  the  cutoff  point  used  (Table  1).  Thus,  differing 
scoring  strategies  do  not  account  for  the  variable  claims  of 
FoxP3  prognostic  significance. 

TIL  can  reside  in  tumor  epithelium,  stroma,  or  both,  and 
this  may  influence  their  prognostic  significance.  Among  the 
58  reviewed  studies,  15  did  not  discriminate  between  the 
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Prognostic  Significance  of  FoxP3+  T  Cells  in  Human  Carcinoma 


Table  1.  Characteristics  of  58  studies 

Poor  Neutral 

FoxP3  prognosis  claim  23  23 

Good 

12 

Specific  antibody 

Clone  42 

1 

Custom 

1 

BioLegend 

3 

1 

Abeam 

6 

4 

1 

mAbcam22509 

2 

1 

1 

Novus  Biologicals 

1 

FJK-16S 

1 

206D 

1 

eBioscience 

2 

236A/E7 

10 

8 

5 

mAbcam2251 0 

2 

1 

PCH101 

2 

1 

259D 

1 

eBio7979 

1 

221D/D3 

1 

Scoring  strategy 

Cutoff  point 

Absence  and/or 

1 

3 

3 

presence 

Mean  cutoff 

1 

3 

1 

Median  cutoff 

16 

11 

5 

Other  cutoff 

5 

6 

3 

Counting  location 

General  count 

5 

5 

5 

Tumor  only 

8 

7 

4 

Tumor  and  stroma 

10 

11 

3 

Tissue  used 

Whole  sections 

16 

17 

7 

Tissue  microarray 

7 

6 

5 

Counting  strategy 

Investigator(s) 

16 

14 

8 

Computer  program 

4 

2 

Not  reported 

7 

5 

2 

Multivariate  correction 

for  stage  or  grade 

Yes:  42 

20 

11 

11 

No:  16 

3 

12 

1 

Multivariate  correction 

for  other  TIL  subsets 

Yes:  8 

2 

3 

3 

No:  50 

21 

20 

9 

Use  of  multiple  markers 

Yes:  8 

4 

4 

No:  50 

19 

19 

12 

Tumor  site 

Hepatocellular 

5 

Cervical 

2 

Head  and  neck 

1 

Pancreatic 

1 

(Continued  on  the  following  page) 

Table  1.  Characteristics  of  58  studies  (Cont'd) 


FoxP3  prognosis  claim 

Poor 

23 

Neutral 

23 

Good 

12 

Renal 

1 

1 

Lung 

1 

1 

Endometrial 

1 

2 

Melanoma 

2 

2 

Breast 

5 

1 

2 

Mesothelioma 

1 

Vulvar 

1 

Oral 

1 

3 

1 

Gastric 

2 

1 

2 

Ovarian 

1 

4 

2 

Bladder 

1 

Colorectal 

6 

4 

NOTE:  Study  N  mean  (range):  219  (30-1,445). 


epithelial  and  stromal  location  of  FoxP3+  T  cells  and, 
instead,  provided  a  general  count;  19  counted  only  FoxP3+ 
T  cells  in  the  epithelium;  and  24  counted  FoxP3+  T  cells 
from  epithelial  and  stromal  compartments  independently 
(Table  1).  Regardless  of  the  location  of  enumerated  FoxP3+ 
T  cells,  a  fairly  even  distribution  was  seen  among  poor, 
neutral,  and  good  prognostic  claims. 

We  next  examined  the  use  of  tumor  tissue  microarrays 
(TMA)  compared  with  whole  sections  (Table  1).  TMAs  were 
used  in  18  of  the  58  studies,  and  prognostic  claims  ranged 
from  poor  (n  =  7),  to  neutral  (n  =  6),  to  good  (n  =  5).  A 
similar  range  of  prognostic  claims  was  seen  with  studies 
using  whole  tissue  sections.  Regarding  cell  counting,  38 
studies  used  manual  counting  by  one  or  more  investigators, 
6  studies  used  a  computer-based  quantification  method, 
and  14  studies  did  not  state  the  counting  method  (Table  1). 
Studies  that  used  manual  counting  showed  an  unbiased 
spread  among  poor  (n  =  16),  neutral  (n  =  14),  and  good 
(n  =  8)  prognostic  claims,  regardless  of  the  number  of 
investigators  who  carried  out  cell  counting.  Definitive  con¬ 
clusions  could  not  be  drawn  regarding  the  use  of  comput¬ 
erized  counting,  as  only  6  studies  used  such  methods,  3  of 
which  involved  colorectal  cancer  (see  below). 

Multivariate  Correction  for  Stage  or  Grade  of 
Disease 

In  principle,  the  density  of  FoxP3+ T  cells  could  reflect  the 
stage  and/or  grade  of  disease,  which  could  influence  prog¬ 
nosis.  Of  the  45  studies  that  correlated  FoxP3+  T  cells  to 
stage  and/or  grade,  25  found  a  significant  association 
between  FoxP3+  T  cells  and  the  stage  and/or  grade  of 
disease,  with  11  reporting  a  P-value  <  0.001  (Supplemen¬ 
tary  Table  SI).  A  potential  confounding  effect  is  that  the 
quantity  of  TIL  can  influence  nodal  staging,  especially  in 
colorectal  cancer  (74).  Nonetheless,  these  studies  support 
the  possibility  that  FoxP3+  T  cells  could  simply  serve  as  a 
marker  of  more  advanced  disease. 
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This  issue  was  addressed  in  42  studies  by  use  of  multi¬ 
variate  models  that  included  stage,  grade,  and  other  clini- 
copathologic  features  (Table  1).  Among  these  studies,  the 
prognostic  significance  of  FoxP3+  T  cells  ranged  from 
poor  (n  =  20),  to  neutral  (n  =  11),  to  good  (n  =  11). 
Notably,  in  4  studies,  FoxP3+  T  cells  were  a  significant 
univariate  prognostic  indicator,  only  to  be  removed  dur¬ 
ing  multivariate  analysis.  Of  the  16  studies  that  did  not 
use  multivariate  analysis,  the  potentially  confounding 
effects  of  stage  and  grade  were  mitigated  in  most  by  the 
fact  that  (i)  FoxP3+  T  cells  showed  no  prognostic  signif¬ 
icance  even  in  univariate  analysis  or  (ii)  only  specific 
stages  or  grades  of  disease  were  included  in  the  study.  In 
summary,  even  though  FoxP3+  T  cells  are  frequently 
associated  with  the  stage  and/or  grade  of  disease,  we 
found  that  this  factor  was  well  controlled  in  most  studies 
and  does  not  account  for  the  variability  of  FoxP3 
prognostication. 

Multivariate  Correction  for  Other  Tumor- 
Infiltrating  Lymphocyte  Subsets 

FoxP3+  T  cells  are  usually  found  together  with  other  TIL 
subsets,  which  can  make  it  difficult  to  discern  their  inde¬ 
pendent  prognostic  effect.  Although  multivariate  analysis 
can  solve  this  problem,  it  requires  that  all  TIL  subtypes 
significant  in  univariate  analysis  be  included  in  the  multi¬ 
variate  model.  In  the  8  studies  that  included  all  TIL  subsets 
in  multivariate  analysis,  the  prognostic  value  of  FoxP3+  T 
cells  ranged  from  poor  (n  =  2),  to  neutral  (n  =  3),  to  good 
(n  =  3;  Table  1).  Thus,  although  the  number  of  studies  is 
low,  it  seems  that  the  prognostic  significance  of  FoxP3+  T 
cells  is  not  solely  attributable  to  the  presence  of  other  TIL 
subpopulations. 

Several  studies  made  prognostic  claims  on  the  basis  of 
the  ratio  of  FoxP3+  T  cells  to  other  lymphocyte  subsets, 
including  CD3+/FoxP3+  (n  =  2),  CD4+/CD25+FoxP3+ 
(n  =  1),  CD68+/FoxP3+  (n  =  2),  CD8+/FoxP3+  or 
FoxP3+/CD8+  ( n  =  18),  CD8+/CCR4+FoxP3+  (n  = 
1),  FoxP3+/CD4+  (n  =  2),  FoxP3+/CD3+/CD45RO+ 
(n  =  1),  and  Granzyme-B+/FoxP3+  [n  =  1;  Supplemen¬ 
tary  Table  SI).  Among  these  28  studies,  prognostic  claims 
for  FoxP3+  TIL  ranged  from  poor  (n  =  12),  to  neutral 
(n  =  11),  to  good  (n  =  5).  Thus,  the  use  of  lymphocyte 
ratios  has  been  inconsistently  applied  and  yielded  incon¬ 
sistent  prognostic  claims. 

Clinical  Significance  and  Publication  Bias 

We  next  evaluated  whether  the  magnitude  of  the  prog¬ 
nostic  effect  was  similar  for  studies  claiming  good  compared 
with  poor  prognosis.  Of  the  58  studies,  32  reported  mul¬ 
tivariate  hazard  ratios  for  overall  survival.  A  funnel  plot 
revealed  no  significant  difference  between  the  magnitude  of 
hazard  ratios  for  studies  claiming  poor  compared  with  good 
prognosis  (Fig.  1).  Furthermore,  there  was  no  evidence  of 
publication  bias,  as  the  studies  were  evenly  distributed 
throughout  the  plot. 


Use  of  Multiple  Markers  to  Define  FoxP3+  T  Cells 

Although  FoxP3  was  originally  thought  to  uniquely 
define  conventional  CD 4+  Tregs  (75),  more  recent  studies 
indicate  that,  in  some  circumstances,  FoxP3  can  also  be 
expressed  by  effector  T  cells  (16,  18).  We  assessed  whether 
studies  that  subdivided  FoxP3+  T  cells  using  a  second 
marker  yield  more  consistent  prognostic  results.  Of  the 
58  reviewed  studies,  50  used  FoxP3  as  a  sole  marker,  which 
resulted  in  variable  prognostic  claims  ranging  from  poor  (n 
=  19),  to  neutral  (n  =  19),  to  good  (n  =  12;  Table  1).  The 
remaining  8  studies  measured  at  least  one  marker  in  addi¬ 
tion  to  FoxP3,  including  CD4,  CD 8,  CD25,  and  C-C  che- 
mokine  receptor  4  (CCR4).  Four  of  these  8  studies  showed 
that  FoxP3+  T  cells  that  coexpressed  a  second  marker  were 
associated  with  poor  prognosis.  The  remaining  4  claimed 
that  the  identified  subset  did  not  have  any  prognostic 
significance.  Of  note,  none  of  the  8  studies  claimed  an 
association  with  good  prognosis. 

On  the  basis  of  the  above  findings,  we  investigated  more 
closely  which  markers  were  used  in  addition  to  FoxP3.  Shah 
and  colleagues  used  2-color  IHC  to  identify  both 
CD4+FoxP3+  and  CD8+FoxP3+  T  cells  in  cervical  cancer. 
Intriguingly,  they  found  CD8+FoxP3+  T  cells  at  a  mean 
number  of  3.32  per  high-power  field  and  CD4+FoxP3+  T 
cells  at  a  mean  number  of  1 1.45  per  high-power  field  (30). 
Thus,  had  they  used  FoxP3  as  a  single  marker,  only  ~75%  of 
the  cells  they  measured  would  have  been  CD4+  T  cells, 
which  underscores  the  fact  that  not  all  FoxP3+  T  cells  are 
conventional  Tregs.  In  another  study,  Watanabe  and  col¬ 
leagues  used  coexpression  of  CCR4  to  delineate  a  subset  of 
FoxP3 -expressing  T  cells  in  oral  cancer  (62).  An  average  of 
58%  of  FoxP3+  cells  were  found  to  coexpress  CCR4. 
Whereas  total  FoxP3+  T  cells  had  no  prognostic  value 
(similar  to  3  other  studies  of  oral  cancer;  refs.  60, 
61,  63),  CCR4+FoxP3+  T  cells  showed  a  highly  significant 
association  with  survival.  These  studies  highlight  the  impor¬ 
tance  of  using  additional  markers  to  account  for  the  het¬ 
erogeneity  of  FoxP3+  T  cells. 

Tumor  Site  and  Subtype 

It  is  conceivable  that  the  biologic  and  prognostic  effect  of 
FoxP3+  T  cells  depends  on  microenvironmental  context,  in 
which  case,  tumor  site  and  histologic  and/or  molecular 
subtype  may  be  important  factors.  Indeed,  when  tumor 
site  was  taken  into  consideration,  we  found  clear  prognostic 
associations  in  some  cases.  For  example,  the  5  studies  of 
hepatocellular  cancer  unanimously  concluded  that  FoxP3+ 
T  cells  are  associated  with  a  poor  prognosis  (Table  1). 
Conversely,  4  out  of  10  studies  investigating  colorectal 
cancer  concluded  that  FoxP3+  T  cells  correlated  with  a 
good  prognosis,  whereas  the  remaining  6  studies  found  no 
prognostic  association.  In  considering  colorectal  cancer, 
Ladoire  and  colleagues  recently  hypothesized  that  the  favor¬ 
able  prognostic  effect  of  FoxP3+  T  cells  may  reflect  their 
ability  to  suppress  tumor-promoting  inflammatory 
responses  to  gut  microbes  (76). 
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Figure  1.  Analysis  of  clinical 
significance  and  publication  bias. 
The  figure  shows  a  funnel  plot  of 
log-transformed  hazard  ratios  (HR) 
compared  with  standard  error  for 
the  reviewed  studies.  Each  symbol 
represents  one  study:  ▼,  poor 
prognostic  claim;  A,  good 
prognostic  claim;  □,  neutral 
prognostic  claim.  Bars  represent 
95%  confidence  intervals. 


In  contrast  to  the  above  examples,  the  prognostic  signif¬ 
icance  of  FoxP3+  T  cells  remains  controversial  in  several 
other  cancers.  In  breast  cancer,  the  reported  prognostic 
effect  of  FoxP3+  T  cells  ranges  from  poor  (n  =  5),  to  neutral 
(n  =  1),  to  good  (n  =  2).  Although  ovarian  cancer  was  one  of 
the  first  tumor  sites  in  which  CD4+  Tregs  were  associated 
with  poor  prognosis  (5),  subsequent  studies  using  FoxP3  as 
a  marker  are  split  among  poor  (n  =  1),  neutral  (n  =  4),  and 
good  (n  =  2)  prognostic  claims.  Similarly,  studies  looking 
at  gastric  cancers  show  a  split  among  poor  (n  =  2),  neutral 
(n  =  1),  and  good  (n  =  2)  prognostic  claims.  For  the 
remaining  10  tumor  sites,  the  number  of  published  studies 
is  insufficient  to  make  definitive  conclusions  about  the 
prognostic  significance  of  FoxP3+  T  cells. 

In  addition  to  tissue  of  origin,  tumors  can  be  classified 
based  on  their  molecular  features,  as  discussed  recently  by 
Ogino  and  colleagues  (77).  Hence,  it  is  conceivable  that  the 
variability  of  FoxP3+  T  cell  prognostication  could  be  attrib¬ 
utable  to  the  inherent  molecular  heterogeneity  within 
tumor  types.  In  support  of  this  idea,  the  prognostic  value 
of  FoxP3+  T  cells  is  stronger  in  mismatch  repair-proficient 
colorectal  cancer  compared  with  mismatch  repair-deficient 
colorectal  cancer  (31).  Similarly,  FoxP3+  T  cells  are  prog- 


nostically  significant  in  estrogen  receptor  (ER)+  but  not 
ER—  breast  cancer  (22,  27).  In  uveal  melanoma,  FoxP3+  T 
cells  provide  prognostic  significance  in  cyclooxygenase-2 + 
cases  (58).  Although  few  in  number,  these  studies  suggest 
that  the  molecular  subtype  of  tumors  may  influence  the 
prognostic  value  of  FoxP3  T  cells. 

Conclusions 

Having  critically  reviewed  the  literature  on  the  prognostic 
value  of  FoxP3+  T  cells,  we  can  make  several  recommenda¬ 
tions  for  future  studies,  (i)  We  recommend  that  prognostic 
marker  studies  follow  a  standard  reporting  structure,  such  as 
the  REMARK  criteria  (78).  (ii)  In  many  cancers,  FoxP3+  T 
cells  are  highly  correlated  with  the  stage  and  grade  of 
disease;  therefore,  it  is  important  to  correct  for  these  and 
other  appropriate  clinicopathologic  factors,  (iii)  FoxP3+  T 
cells  are  invariably  found  with  other  lymphocytes;  there¬ 
fore,  all  TIL  subsets  with  prognostic  value  should  be  includ¬ 
ed  in  multivariate  models,  (iv)  The  use  of  multiple  markers 
to  identify  functional  subsets  of  FoxP3+  T  cells  can  lead  to 
greater  clarity  about  their  prognostic  value,  (v)  The  prog¬ 
nostic  value  of  FoxP3+ T  cells  seems  to  depend  significantly 
on  tumor  site  and  possibly  molecular  subtype,  suggesting 
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that  the  biologic  properties  of  FoxP3+  T  cells  are  influenced 
by  the  tumor  microenvironment  in  which  they  reside. 
Overall,  this  study  provides  a  cautionary  note  for  the  con¬ 
cept  of  depleting  FoxP3+  cells  from  patients  with  cancer  as  a 
means  to  enhance  tumor  immunity.  Our  findings  suggest 
that  this  strategy  may  be  beneficial  for  some  tumor  sites 
(e.g.,  liver)  but  detrimental  to  others  (e.g.,  colorectal). 
Improved  understanding  of  the  different  FoxP3+  T  cell 
subsets  in  human  cancer  will  likely  enable  the  development 
of  more  precise  and  effective  immunotherapies. 
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